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CHAPTER  I 


Introduction 


The  diffraction  of  scalar  (acoustical)  waves  by  the  tip  of  an  eHiptical  perfectly 
conducting  cone  was  studied  by  Krauss  and  Levine  [1].  Satterwkite  and  Kouy- 
oumjian  [2]  examined  the  vector  electromagnetic  problem  and  presented  a  Green’s 
dyadic  for  a  source  radiating  in  the  pr>  ^ence  of  an  angular  sector.  However,  their 
solution,  expressed  in  terms  of  non^closed  form  Lame  functions,  is  cumbersome 
for  numerical  calculations.  Furthermore,  so  far  it  has  not  appeared  possible  to 
asymptotically  identify  a  “corner  diffraction  coefficient”  from  this  eigenfunction 
representation. 

Recently,  Burnside  and  Pathak  |3]  proposed  a  corner  diffraction  coefficient 
which  succesfuliy  predicted  the  corner  effect  of  numerous  plate  structures.  Their 
solution  is  based  on  the  asymptotic  evaluation  of  the  radiation  integral  involving 
the  equivalent  currents  that  would  exist  in  the  absence  of  the  corner.  A  corner 
diffraction  term  is  then  established  by  heuristically,  but  at  present  empirically, 
modifying  the  final  result.  Sikta  (4]  modified  the  spread  factor  of  the  diffracted 
field  in  [3]  and  applied  a  limiting  process  to  derive  the  wave  diffracted  by  the 
corner  and  propagating  along  one  of  the  edges  of  a  plane  right  angular  sector. 
By  introducing  an  empirically  established  “reflection  coefficient”  he  utilized  his 
edge  wave  corner  diffraction  coefficient  in  the  calculation  of  the  double  and  triple 
diffraction  by  two  adjacent  corners  of  a  flat  plate  structure. 
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The  equivalent  current  approach  is  adopted  herein  to  study  the  vertex  diffrac¬ 
tion  of  an  ele  troinagnetic  wave  guided  along  one  of  the  edges  of  a  semi-infinite 
wedge.  The  geonletry  of  the  problem  is  illustrated  in  Fig.  1.  The  dipole  source  is 

t 

rs  diating  in  the  close  vicinity  of  the  edge  and  excites  a  paraxial  field  guided  by  one 
o'/  the  edges  of  the  trihedron  -  hence  the  term  “edge  wave".  Explicit  expressions  of 
'  this  field  are  derived  in  Chapter  II  based  on  the  limiting  behavior  of  the  Green’s 
dyadic  for  an  infinite  wedge  [7]  and  small  distances  of  the  dipole  from  the  edge,  and 
the  canonical  solution  (Fresnel  integral  representation)  to  the  half  plane  problem 
for  plane  wave  grazing  incidence.  A  simple  edge  diffraction  coefficient  can  then  be 
established  valid  for  sinall  distances  of  the  source  from  the  edge  and  for  paraxial 
field  calculations. 


Figure  1:  An  infinitesimal  electric  dipole  source  radiating  in  the  close  vicinity  of 
one  of  the  edges  of  a  semi-infinite  wedge. 

In  Chapter  III  a  Physical  Optics  approximation  of  the  vertex  diffracted  field 
is  presented,  based  on  the  currents  that  would  be  induced  by  the  edge  guided 
wave  for  an  infinite  wedge  which  is  then  truncated.  The  asymptotic  evaluation 


2 


of  th«  surface  currents  radiation  integral  appropriately  encounters  the  edge  wave 
singularity  consistent  with  Meixner’s  edge  condition  [5]. 

The  concept  of  equivalent  currents  is  generally  examined  in  Chapter  IV.  The 
end  point  contribution  to  the  rar  'ntion  integral  of  these  currents  “excited”  by  the 
impinging  edge  wave  is  interpreted  as  a  vertex  diffraction  term.  Two  types  of 
equivalent  currents  for  arbitrary  aspects  of  observation  is  shown  to  yield  compara¬ 
ble  results.  Edge  wave  vertex  and  edge  wave  edge  diffraction  coefficients  can  then 
be  established  by  an  empirical  modification  of  the  asymptotic  field  expressions. 
The  validity  of  the  approach  is  confirmed  via  comparison  with  moment  method 
results  and  pattern  measurements  for  a  small  dipole  radiating  in  the  close  vicinity 
of  one  of  the  edges  of  a  polygonal  plate. 

Every  structure  considered  in  the  present  work  is  assumed  perfectly  conduct¬ 
ing.  Furthermore,  an  exp  (ju>of )  time  dependence  has  been  adopted  and  suppressed 
in  the  following  analysis. 
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CHAPTER  II 


Edge  Waves 


2.1  Introduction. 

The  term  “edge  waves”  in  the  present  work  defines  waves  propagating  along 
the  edge  of  a  (possibly  curved)  wedge.  The  edge  wave  is  actually  a  form  of  a 
maxwellian  field  guided  by  the  edge,  and  exhibiting  the  proper  singularity  in  ac¬ 
cordance  with  Meixner’s  edge  condition. 

As  shown  later  in  this  chapter,  such  a  paraxial  singular  field  can  be  excited 
either  by  a  plane  wave  at  grsizing  incidence,  or  by  a  dipole  radiating  in  the  close 
vicinity  of  the  edge.  The  vertex  of  a  terminated  edge  illuminated  by  a  plane  or 
spherical  wave  can  also  excite  an  edge  wave.  Independently  of  the  excitation, 
but  sufficiently  far  from  its  source,  the  edge  wave  behaves  and  can  be  treated 
as  a  ray  optical  field.  However,  application  of  ray  optical  techniques  (UTD)  is 
not  straightforward,  mainly  due  to  the  singular  behavior  of  the  paraxial  fields. 
An  asymptotic  high  frequency  approximation  which  encounters  this  peculiarity  is 
attempted  in  Chapters  III  and  IV. 

Two  types  of  excitation  of  an  edge  wave  are  examined  in  this  chapter;  a  plane 
wave  grazing  the  edge  of  a  half  plane  and  an  elementary  dipole  radiating  in  the 
near  vicinity  of  the  edge  of  an  infinite  wedge.  The  limiting  behavior  of  the  fields 
predicted  by  the  canonical  solution  for  the  half  plane  problem  and  the  leading  term 
of  the  power  series  expansion  of  the  Green’s  dyadic  with  respect  to  the  numerical 
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distance  k^p'  of  the  point  source  from  the  edge  (where  ko  =  2n/\  denotes  the 
free  space  wavenumber)  are  examined.  The  two  cases  are  related  via  reciprocity 
so  that  a  simple  edge  diffraction  dyadic  coefficient  can  be  readily  established  and 
applied  in  the  approximation  of  the  patterns  of  complex  sources  radiating  in  the 
close  vicinity  of  the  edge  or  in  paraxial  field  calculations. 

2.2  The  canonical  solution  to  the  half  plane  problem.  Some  limiting 
cases. 


A  perfectly  conducting  half  plane  is  illuminated  by  an  arbitrarily  polarized 
plane  wave,  as  depicted  in  Fig.  2.  For  convenience,  the  origin  of  our  reference  frame 
has  been  chosen  to  coincide  with  the  origin  Q*  of  the  Keller  cone  of  diffracted  rays, 
defined  by  the  angle  of  incidence  and  the  observation  point  2).  Since 

the  only  dependence  of  the  total  field  on  the  z  coordinate  is  incorporated  into  the 
factor  exp{— jk^z  cos  the  problem  is  essentially  quasi  two  dimensional.  It  can 
be  easily  shown  that  the  p  and  (^directed  field  components  can  be  expressed  in 
terms  of  the  Ez  and  Hz  components  via  the  equations 


V1E2  ,  zxVtHz 

Et  -  jkz  ,  2  _  l2  .2  _  ,2 

Wfl  "-z  *'0 

ur  I  ^tHz  2  X  ^tEz 

Ht  =  -jkz  7-5 - -Z  -  JWQCO 


k^  -  fc2 
(Vq  ">2 


kl-kj 


(2.1a) 

(2.1b) 


where  kz  =  fco  ^osfi\  wq  =  kaf  y/pQtQ  and  pQ,  eg  are  the  permeability  and  the 
permittivity  of  the  free  space.  The  operator  V(  is  defined  by 


—  A  ,  d  -  1  5 


(2.2) 


Sommerfeld’s  canonical  solution  to  the  half  plane  problem  can  be  expressed 
in  a  compact  form  in  terms  of  the  modified  Fresnel  Integral.  In  particular,  we  have 
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Figure  2:  Half  plane  illuminated  by  an  arbitrarily  polarised  plane  wave. 


F, 

Ei(Q^) 

[  «i(Qn 

^  exp(-jto*) 


(2.3) 


where  we  substituted 


=  A'_  ( -  y/ 2kop  sin  /?'  cos  ) 


(2.4) 


and  8  can  be  expressed  in  terms  of  the  coordinates  (p,^,z)  of  the  observation  point 
R  as 


s  =5  p  sin  +  z  cos  0' 


(2.5) 


The  function  A'-(x)  is  the  modified  Fresnel  Integral  defined  by  [6] 

A_(t)  =  F_(i)  exp{>(x^  +  ir/4)}  (2.6) 

y/ir 

F_(x)  =  J  txp(-jf^)dt  (2.7) 
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Incorporating  eqs.  (2.3)-(2.7)  into  eqs.  (2.1a)-(2.1b)  and  performing  the  dif¬ 
ferentiations  one  can  readily  obtsdn  a  complete  representatior  of  the  transverse  (to 
the  z  direction)  components  of  the  electric  and  magnetic  fields.  Namely,  one  finds 


cos/3'  [u‘  sin(<^  -  <f>)  -  u*"  sin{<i>  +  <^')] 

+E^^,{Q^)  [u‘  coa(<^  -  <f/)  -f  u’’  cos{(f>  +  4>')] 

,  2exp{-jiT/4)  4> 

-I —  cos  — 

\j2irkQp  sin  ^ 

•{i5^(Q®)  cos/3'sin^  -I-  cos  y]|  cxp{-jkos) 

(2.8b) 

while  the  transverse  magnetic  field  Hf  can  be  obtained  by  using  duality. 

Of  particular  inter'^st  in  our  work  is  the  limiting  behavior  of  the  above  expres¬ 
sions  for  small  values  of  the  parameter 

e(p;/3')  =  fcop  sin/?'  (2.9) 

Eqs.  (2.8a)-(2.8b)  can  then  be  represented  by  a  power  series  expansion  of  the  form 

E(p,4:,z)=  +  (2-10) 

V*  i=0 

As  ({p\0')  0  one  may  write 
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E,(p,^,z)  =  0(<®) 


(2.11) 


=  \/ nkohh^0‘  ^  ^  ‘ 

•  l^^(P')  co«/^'  «n  J  +  -S^J'CQ*)  cos  y  exp(-;fco«) 
+d(c®)  (2.12) 

whfr«  d(e^)  denotes  the  remainder  in  the  power  series  expansion  and  its  order  with 
respect  to  e. 


Figure  3:  An  infinitesimal  electric  dipole  source  radiating  in  the  close  vicinity  of 
the  edge  of  an  infinite  perfectly  conducting  wedge. 

At  grazing  incidence  the  leading  term  in  the  power  series  expansion  of 
r)  represents  a  wave  guided  along  the  edge  of  the  half  plane  (an  edg'- 
wave)  g’ven  by 

E'“'(p,^,r)  =  exp(-Jw/4M/’ »i»^  +  ^cos^) 
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(2.13) 


where 

cy  .  •(4«?')««(^72)  +  £;i((?')co8(^72)]  .  ^ 

S(p,4>)  = - 7^-“ -  (2.14) 

y/tinp' 

is  a  source  related  factor.  The  wave  described  by  eqs.  (2.13),  (2.14)  satisfies  the 
wave  equation,  the  Meixner’s  edge  condition  and  dominates  the  field  in  the  paraxial 
region. 

2.3  The  limiting  behavior  of  the  dyadic  Greenes  function  for  a  dipole 
radiating  in  the  close  vicinity  of  the  infinite  edge  of  a  wedge. 

Let  us  consider  the  more  general  reciprocal  problem  of  a  point  source  = 
j?#(f  —  r')  radiating  in  the  presence  of  a  perfectly  conducting  infinite  wedge  as 
depicted  in  Fig.  3.  The  field  produced  by  this  dipole-wedge  configuration  can  be 
formally  written  as 

£(f)  =  r' ) .  p«  (2.15) 

where  F(r,r')  is  the  Green’s  dyadic  for  the  wedge.  An  eigenfunction  series  repre¬ 
sentation  of  r(f,r’)  is  pven  in  |7). 

For  convenience  the  source  point  is  located  at  5(p\^\0) ,  (r’  =s  p\P'  =  ir/2). 
Furthermore,  let  kor  >>  1,  so  that  the  spherical  Hankel  functions  of  the  second 
kind  involved  in  the  series  representation  assume  their  asymptotic  form.  The 
limiting  form  of  the  Green's  dyadic  for  small  values  of  the  parameter  e  =  kop'  sin^ 
is  of  interest  here.  After  some  tedious  manipulation  it  can  he  shown  that 

f  (f, r' )  =  C{u)  {kop  sin  p' )"" *  (d  cos p sin  ^  cos  u4>) 

{p'  sini/^’-f  cosi/^7 — ^  '  +0(«®)  (2.16) 

T 

with  the  constant  factor  C{v)  defined  by 
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(2.17) 


and 


C(v\  =  d-  2)  gxp{i(i/  +  l)ir/2)} 

'  ’  x/923*'+*r(i/)r(i/  +  i)r(i/  +  3/2) 


1  2n  - 

1/  =  -  ;  n  =  - 

n  w 


(2.18) 


where  (WA)  denotes  the  wedge  angle.  Hence,  the  total  field  for  a  general  wedge  is 
written  as 


£(f)  =  ju)ofioC(t')(kop'Bin/J)*'~^(0  cos/?  sint/^  +  ^cosi/^) 

It'  i!  t  t  ^^P(~J^0^)  I 

•(p^  sin  +  p^f  cos  1/0  ) - +  C/(e  ) 


(2.19) 


which,  in  the  case  of  the  half  plane,  reduces  to  the  following  expression: 

VtoZoexp(>»/4)  I 

= - 

/?  cos^8in(0/2)  4  0  co8(0/2)  exp(-;feor) 

\/sin^  r 


(2.20) 


where 


>4'(p',^') 


p',  sin(072)  +  pj*  cos(072) 


(2.21) 


is  a  constant  source  factor.  The  remainder  0(e®)  in  eq.  (2.19)  can  be  expressed  in 
a  closed  form  only  for  the  case  of  the  half  plane.  Specifically,  it  can  be  obtained 
via  reciprocity  from  eqs.  (2.8a)-(2.8b)  after  subtracting  the  edge  wave  term.  In 
the  case  of  a  general  wedge  angle  the  remainder  is  essentially  a  power  series  of  the 
small  parameter  c,  the  coefficients  of  which  can  be  derived  from  the  eigenfuction 
series  representation  of  the  Green’s  dyadic. 

The  lower  order  term  in  the  RHS  of  eq.  (2.19)  dominates  in  the  paraxisi 
region  (edge  wave)  and  reveals  the  strong  couplinx  between  the  dipole  source  and 
the  edge.  It  is  a  spherical  ray  optical  wav^  that  satisfies  Maxwell's  equations. 
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Cltarly,  the  behavior  of  the  electric  field  intensity  of  the  edge  wave  is  independent 
of  the  orientation  of  the  source. 

2.4  The  half  plane  case.  : 


As  pointed  out  in  the  previous  section  an  exact  closed  form  representation  of 
the  Fraunhofer  region  field  produced  by  the  dipole  ^  radiating  in  the  presence 
of  an  infinite  half  plane  can  be  obtained  by  applying  the  reciprocity  theorem  to 
eqs.  (2.8a)-(2.8b).  For  small  distances  of  the  dipole  from  the  edge  one  can  retain 
only  the  two  leading  terms  in  the  power  series  expansion  of  the  modified  Fresnel 
integrals  involved  in  the  canonical  solution.  It  turns  out  that  the  total  radiated 
field  ran  be  described  by  the  superposition  of  three  terms: 

i)  The  incident  field  produced  by  the  dipole  in  the  absence  of  the  half  plane 
multiplied  by  a  factor  of  1/2, 

ii)  The  field  produced  by  the  image  of  the  dipole  source  with  respect  to  an 
infinite  plane  multiplied  by  a  factor  of  1/2, 

iii)  An  edge  associated  field,  which  may  be  treated  as  an  edge  diffracted  wave 
(emanating  from  Q')  in  a  UTD  sense. 

It  should  be  noted  that  the  components  i),  ii)  illuminate  the  whole  space  and 
not  only  the  regions  bounded  by  the  Geometrical  Optics  incident  and  reflected 
shadow  boundaries.  Under  this  consideration  the  total  field  is  uniform  everywhere 
and  satisfies  the  appropriate  boundary  conditions  on  the  conducting  half  plane.  In 
a  matrix  notation  one  can  write  for  the  radiated  field 


I  Ep{T) 

1 

;  £}■•(?) 

jkoZo  cxpl-j  tro’') 


8x 


[exp{;ibo^'sinAfcos(^  ~ 
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—  co8^co8(^  —  —  co8^8in(^  —  lin/9 

iiu(^  -  4>')  -  co8(^  -  ^')  0 

+  exp{  j  kop'  sill  /9  co8(^  + 

co8/9co8(^  +  ^')  -cos^iin(^  +  ^')  —  lin/J 
-  •in(^  +  4>')  -  coi(^  +  ^')  0 


(2.22) 


in  which  E**^'(f)  denotes  the  edge  associated  term.  The  latter  can  be  readily 
obt«ined  from  eqs.  (2.13),  (2.14)  or  eqs.  (2.20),  (2.21).  By  empirically  generalizing 
those  results  the  diffracted  field  can  be  expressed  in  a  more  familiar  notation  by 
the  following  equation: 


2exp(-jir/4) 
y/2irko  sin  0o 


tin  ^  sin  |  cos  ^  f 
sec  0Q  sin  ^  cos  ^  cos  ^  cos  ^ 


(2.23) 


The  various  parameters  involved  in  the  above  expression  are  shown  in  Fig.  3.  Also, 
E'Z  .i(Q')  nre  the  ray  optical  components  of  the  total  field  incident  at  the  point 
and  produced  by  the  dipole  source,  namely 


1  I  r^' 


(2.24) 
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with  /3q  =  /?o  denoting  the  diffraction  angle.  It  should  be  mentioned  that  eq.  (2.23) 
reduces  to  the  correct  limiting  expression  of  the  field  in  the  paraxisd  region  of  the 
edge  {I3q  — »  0),  exhibiting  the  appropriate  edge  singularity. 

2.5  Discussion  and  numerical  results. 

The  accuracy  of  the  “edge  wave”  solution  given  by  eqs.  (2.23)-(2.24)  in  pre¬ 
dicting  the  field  for  small  distances  of  a  point  dipole  source  from  the  edge  of 
a  half  plane  is  illustrated  in  Figs.  4-9,  in  comparison  with  the  exact  Sommer- 
feld  (Fresnel  integral)  representation,  for  several  azimuthal  angles  <)>  of  the  ob¬ 
servation  point  which,  for  simplicity,  is  removed  at  infinity.  The  infinitesimal 
source  has  a  dipole  moment  radially  directed  and  it  is  located  at  the  points 
(p\  <!,•)  =  (0.025Ao,  135“),  (0.05Ao,  75®),  (0.075Ao,  180®),  in  terms  of  the  asso¬ 
ciated  system  of  cylindrical  coordinates  centered  at  the  projection  of  the  point 
source  onto  the  edge  of  the  half  plane.  The  amplitude  pattern  of  the  /?—  and 
(^—directed  field  components  is  then  plotted  against  the  elevation  angle  of  the  ob¬ 
servation  point.  Comparison  with  the  exact  representation  of  the  total  field  shows 
that  a  distance  p'  <  O.OSAq  is  adequate  for  the  total  field  to  be  accurately  de¬ 
scribed  (within  1.0  dB  of  magnitude)  by  eqs.  (2.23)-(2.24).  As  a  matter  of  fact, 
the  accuracy  is  improved  for  smaller  distances  of  the  dipole  or  the  field  point  (i.e., 
when  grazing  incidence  is  approached)  from  the  edge  of  the  half  plane.  In  the 
latter  case  the  total  field  exhibits  a  singular  behavior  which  is  consistent  with  the 
edge  condition. 

Clearly,  the  analysis  of  the  waves  associated  with  the  edge  of  an  infinite  wedge 
fails  in  the  case  of  finite  or  semi-infinite  edges.  First,  the  edge  wave  term  should  be 
suitably  modified  so  that  the  singularity  in  the  extension  of  the  edge  is  eliminated. 
Second,  the  corner  effects  should  be  incorporated  in  the  total  solution.  Asymp- 


13 


totically,  the  end  point  efTect  is  additively  introduced  into  the  total  solution  and 
corresponds  to  a  corner  diffracted  ray.  On  the  other  hand,  the  edge  wave  singular¬ 
ity  is  compensated  multiplicalively  with  the  introduction  of  the  proper  transition 
dyadic.  Within  this  context,  for  a  dipole  radiating  in  the  close  vicinity  of  a  ter¬ 
minated  edge  and  sufficiently  far  from  the  corner  (in  the  angular  sector  case),  the 
total  field  can  be  written  as 

jg*  =  +  £'*"•  f  -I-  (2.25) 

where  is  half  the  field  produced  by  the  dipole  and  its  image  with  respect  to 
the  plane  of  the  plate  in  free  space,  is  the  sum  of  the  corner  and,  possibly,  the 
terminating  edge  diffracted  field  and  7  is  a  transition  dyadic  to  be  determined. 
An  empirical  definition  of  this  transition  dyadic  is  ^ven  in  Chapter  IV. 
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Fisnw  «:  P«  NsiMi  ^.dirwtnl  wan  of  the  ioRnitooinul  Figo»o  9:  <ioW  «*•*»**  ^R'WrtoH  woof  of  tko 

dipole  tource  of  Fig.  3  ot  m  uiiiiulhal  angle  4  =  1«S*  for:  dipole  eourre  of  Fig.  3  at  an  aiimnlhal  angle  d  = 

A)  p'  =  0.025A«,  d'  =  I35*,  B)  p'  =  O.OSAo,  d'  *  75*,  A |  p'  =  0.025Ao.  d'  =  135®.  B)  p  =  0.05A«.  d  *  75  , 

C)  p*  =  0.075A,.  d'  =  laO®.  C)  p'  =  0.075AO.  d'  =  1«»®. 


CHAPTER  III 


A  Physical  Optics  Approximation  of  the  Edge  Wave  Vertex  Diffracted 

Field 


3.1  Introduction. 

One  of  the  simplest  approaches  to  high  frequency  scattering  problems  is  the 
Physical  Optics  approximation.  This  method  assumes  the  surface  of ‘a  perfectly 
conducting  scatterer  is  locally  plane  and  approximates  the  total  field  at  the  points 
of  the  surface  illuminated  by  the  incident  wave  with  the  superposition  of  the  inci¬ 
dent  and  reflected  field.  The  field  is  considered  tero  at  the  shadowed  part  of  the 
surface.  The  scattered  wave  is  then  evaluated  via  the  radiation  integral  of  the  in¬ 
duced  surface  currents.  Normally,  Physical  Optics  fails  to  adequately  approximate 
the  scattered  field  in  observation  directions  sufficiently  away  from  the  specular  and 
backscatter  directions. 

The  study  of  the  edge  wave  corner  diffraction  mechanism  with  the  Physical 
Optics  approach,  however,  is  conceptually  different.  The  excitation  dipole  is  lo¬ 
cated  in  the  close  vicinity  of  the  edge  of  a  finite  wedge  and  sufficiently  far  from 
its  vertex  as  shown  in  Fig.  10.  The  induced  surface  currents  are  approximated  by 
the  currents  that  would  flow  over  the  surface  of  an  infinite  wedge.  The  end  point 
contribution  to  the  radiation  integral  is  then  interpreted  as  a  corner  diffracted 
wave. 

This  approach  is  expected  generally  to  yield  an  acceptable  approximation  of 
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the  diffracted  field  whenever  the  edge  wave  current  flow  lines  are  not  significantly 
distorted  in  the  vicinity  of  the  vertex.  Moment  method  results  in  general  justify 
this  assumption. , 


t 


Figure  10:  Geometry  for  the  edge  wave  edge  and  vertex  diffraction  problem. 

S.2  Formulation  of  the  Phyaical  Optica  aolution. 

For  convenience,  let  us  restrict  our  attention  to  face  1  ("O"  face)  of  the  wedge. 
The  result  for  face  2  face)  can  be  readily  obtained  by  means  of  a  simple 
transformation.  The  field  associated  with  the  currents  ii(^)  flowing  over  the  face 
1  surface  can  be  evaluated  via  the  radiation  integral 

=  (3.1) 

where  $e  i>  the  vector  from  pointing  to  the  observation  point,  R  =  ($(  -  r')/R 
and  R  =  Hsc  -  r'||  (Fig.  11).  The  intecratinn  takes  place  over  the  truncated  face 
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of  the  wedge.  It  is  presumed  that  ji(r*)  can  be  adequately  approximated  by  the 
actual  induced  currents  as  if  the  wedge  was  infinite,  namely 

xJ?'-(P)|s|  (3.2) 


Figure  11:  Geometry  for  the  asymptotic  evaluation  of  the  radiation  integral  of 

the  surface  edge  wave  currents. 

where  H^*'{r')  is  the  magnetic  field  produced  by  the  dipole  radiating  in  the  close 
vicinity  of  the  edge  of  an  infinite  wedge.  The  end  point  (vertex)  contribution  to 
the  surface  radiation  integral  of  eq.  (3.1 )  is  of  interest  here.  For  convenience,  let  us 
introduce  the  oblique  system  of  coordinates  associated  with  the  terminating  edge 
(edge  (1)  in  Fig.  10),  defined  by  the  unit  vectors 

V  =:  z  ;  #1  =  i  sinoi  -  i  cosoi 

so  that; 

rs(|sino]  ;  r  =  tt-<icosaj 
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In  the  above  equations  Q|  denotes  the  angle  formed  by  the  guiding  edge  and 
the  terinJnatiiig  edge  (edge  (1))  of  face  1  of  the  wedge.  The  Jacobian  determinant 
of  the  above  transformation  is  then  found  to  be  equid  to 


a(tt,fi)| 


sinai 


and  eq.  (3.1)  results  into  the  following  expression: 


(3.3) 

The  first  term  in  the  RHS  of  eq.  (3.3)  expresses  the  guiding  edge  (edge  (0)  of 
Fig.  10)  associated  wave  (edge  wave)  multiplicatively  corrected  by  the  transition 
dyadic  T{sc),  due  to  the  semi-infinite  rather  than  infinite  nature  of  the  structure, 
whereas  the  second  term  incorporates  the  edge  (1)  effect  as  well  as  the  vertex 
contribution.  Thereafter,  our  attention  is  restricted  to  the  second  term  only. 

The  double  integral  in  eq.  (3.3)  can  be  reduced  to  a  line  integral  along  the 
terminating  edge,  by  considering  the  end  point  contribution  with  respect  to  the 
variable  u.  A  non-uniform  asymptotic  approximation  readily  follows  after  integra¬ 
tion  by  parts.  Next,  a  uniform  representation  is  derived  by  heuristically  introducing 
the  suitable  transition  function  to  the  final  result. 

For  large  distances  of  the  dipole  from  the  terminating  edge,  using  the  results 
of  section  2.3,  one  observes  that 


Jl(«3<l) 


in  which 


;  „.o , 

- - ^  f  gjj]  00  4  ;  (-os  1^0  ) 

**0 


(3.4) 
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sin^O  =  fisinai/ro  ,  cos^O  =  (^  ~  +  *e)/*’0* 


At(p\<l>)  =  Binuif)  +  pj,  cosi/^'), 

t 

*{.  is  the  distance  of  the  dipole  source  from  the  vertex,  and  the  constant  C(i/)  is 
defined  by  eq.  (2.17).  It  is  renunded  that  i/  =  l/n,  where  (2  -  n)x  is  the  wedge 
angle.  Also,  -<c)  deternunes  approximately  the  location  of  the  dipole  with 

respect  to  the  corresponding  cylindrical  system  of  coordinates.  Thus,  eq.  (3.3)  may 
be  written  explicitly 


noo  ^  ^ 

^  A  X  A  X  (f  sin  ^0  +  *  cos  ^o)  wn*^"*  ^0 


exp{-jfco(^-»-co)} 

4xi2r0 


dudti 


(3.5) 


whence,  integrating  by  parts  with  respect  to  the  variable  u,  one  obtains  the  (non- 
uniform)  asymptotic  approximation: 


Elin  ^  E\'‘'{sc)  T{f)  -  jk^Zo  mx'' ai 

'  R  X  R  X  \iii  sinoi  -I-  —  <i  cosoi)]fj~* 

*,  (3.6) 

».(«  +  ’-o)l.=o 

It  is  well  known  that  Physical  Optics  does  not  predict  the  correct  edge  diffrac¬ 
ted  field  away  from  the  geometrical  optics  shadow  boundary.  Therefore,  it  would 
be  reasonable  at  this  point  to  be  restricted  into  the  vertex  contribution  of  the 
integral  involved  in  eq.  (3.6).  Note  that  the  latter  is  essentially  a  sum  of  integrnl'^ 


of  the  form: 

^  txp{jkg{1))dt 
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where  G{t)  is  a  slowly  varying  function  of  t  and  ??(p)  <  0.  Stich  types  of  integrals 
are  evaluated  asymptotically  in  Appendix  A  and  result  in  expressions  involving 
the  parabolic  cylinder  function  of  order  p.  Thus,  the  vertex  diffracted  field  associ¬ 
ated  with  face  1  of  the  wedge  and  predicted  by  the  end  point  contribution  to  the 
radiation  integral  of  the  Physical  Optics  currents  can  be  expressed  as  follows: 

.in;- a,  C{u)AUp',4’') 

47racSc'"^^(l -COS/3)  ^ 

— (sin/3cosai  -I-  sinaj  cos/3cos^)/®^_j(A!o)]4 

-|-sinoisin^/®^_l(A:o)^|  •  (3.7) 


In  the  above  equation,  Ip(ko)  expresses  the  end  point  contribution  to  the  integral 

j^p(fco)  =  ^  exp{-jko(R  +  ro)|u=o}  dt 

which  is  examined  in  detail  in  Appendix  A.  For  the  specific  semi-infinite  wedge 
geometry  under  consideration,  the  integral  Ip(ko)  can  be  approximated  by: 


/“(to)  »  r(-p)*5exp(-j>r/2)elp{-ito(»c  +  4)} 

[sinoj  sin/3co8^ -1- cosai(l  —  cos)9)]  ^ 

in  which 


(3.8) 


“1  ~^l)  \/(»c  +  >c)  -  (»i  +  »l) 


(3.9) 


and  Sj,  Ji  denote  the  distance  of  the  dipole  source  point  and  the  receiver  from  the 
origin  (?j  of  the  Keller  cone  of  diffracted  rays  from  the  terminating  edge  (edge  (1 ) 
in  Fig.  10).  The  branch  of  the  bracketed  expression  in  eq.  (3.8)  is  chosen  according 
to 


(-1)  P  =  exp(-jpir) 
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Recall  that  cq.  (3.7)  furnishes  a  non-uniform  approximation  of  the  corner 
diffracted  field  associated  with  face  1  of  the  wedge.  Moreover,  this  becomes  evident 
from  the  singularity  of  the  field  as  the  receiver  approaches  the  z-axis  {0  -*  0).  A 

I 

uniform  approximation  incorporates  the  coupling  between  the  vertex  and  the  point 
source  and  mathematically  implies  the  multiplicative  introduction  of  the  function 

^ 

where  F(-)  is  the  familiar  edge  transition  function  of  UTD.  The  above  expression 
can  be  derived  from  the  uniform  asymptotic  approximation  of  integrals,  the  inte¬ 
grands  of  which  involve  a  phase  function  presenting  a  stationary  phase  point  (in 
our  case  the  dipole  source)  in  the  vicinity  of  an  end  point  (the  vertex).  It  should 
also  be  noted  that  the  above  transition  function  appears  in  the  edge  wave  corner 
diffraction  coefficient  obtained  in  (4). 

3.3  The  Physical  Optics  vertex  diffracted  wave. 


Eq.  (3.7)  along  with  eqs.  (3.8),  (3.9)  can  be  used  to  identify  a  vertex  diffraction 
coefficient  associated  with  the  Physical  Optics  currents  flowing  over  the  wedge 
surface.  Retaining  only  the  dominant  term  in  the  expression  of  the  vertex  diffracted 
wave  (dominant  with  respect  to  the  parameter  2iboS()  one  obtains 

4irsc«e 

•  cot  I  F(2toi,  .in*  0/2)  ( I  ; 

2  ( (cot  oi(l  —  cos/j)  +  sinj 

. 


i  0  COS 

(3.ini 


(cot  021 1  -  cos  0 )  -f  sin  0  cos(  nn 
Eq.  (3.10)  can  be  rewritten  in  a  more  familiar  notation  by  employing  the  follow¬ 
ing  expressions  for  the  incident  fieM  «ni  th'-  produced  by  the  dipole 
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radiating  in  free  $pace: 


I « gLbJfp.»^)  j 

4«?')  I  A, 


Then, 


(3.11a) 


txp(juir/2)T{u)C{u)  -t-  ) cos 

^0  tin*"*'  0' 


•col  2  F{2kQLeKn 


‘”i>{ 


_ ^i/(feooS(tc;<c)l _ 

(cot  02(1  —  COB0)  +  tin  j9  cot(nit  —  ^)]*' 


^i/lfcogl(<c;<c)] 
(cota|(l  —  cot/9)  -f  tin/? cot 

•♦4)1 _ 1  exp(-i*0«c) 


(3.11b) 


wliere  oi ,  02  are  the  angles  formed  by  the  terminating  edges  (1)  and  (2)  and  the 
guiding  edge  (0)  for  faces  1  and  2  respectively  and  the  parameters  01,2  are  defined 


«1.2(^:;4)  =  -01,2“ ^1,2)  |y(tc  +  4)-(4,2  +  'U)l 

with  ^i^2’^1.2p  ®*  before,  denoting  the  distances  between  the  dipole  or  the  receiver 
and  the  origin  of  the  Keller  cone  of  diffracted  rays  from  the  edges  (1)  and  (2)  of 
faces  1  and  2  respectively,  while  0i  ,02  the  elevation  angles  of  the  observation 
point  in  terms  of  the  edge  (1)  and  (2)  fixed  coordinate  systems.  The  transition 
function  fl,,('),  which  it  examined  in  detail  in  Appendix  B,  is  defined  by 

Fl„(|rl)=exp(ii/ir/4)(2|Tl)"/2exp(.HTl/2)r_„{exp(>ir/4)v^)  (3.13) 


2b 
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For  th«  particular  case  of  the  plane  angular  sector  (n  =  2,  i/  =  1/2  and  01  =  02  = 
a,  tti  =  02  =  a)  eq.  (3.11b)  reduces  to 

^  •  1  £yQ^).in(^'/2)  +  Eyty)co.(//2) 

^  y/2itko  ^sin/?' 

•cot|f(2iol.iin*|) 

exp(->to»c)  ,, 

■■  ■■  ■■  ■-  ■  -  (3.18) 

yeot  o  (1  —  cos 0)  +  sin/?  cos 4> 

and  for  the  right  angular  sector  (a  =  ir/2)  it  simplifies  to  the  following 

_ 1  4(Q^)  sin(^72)  -f  £;,(Q^)  cos(^72) 

^  \/2  ir  fco  ^sin 

•  cot  ^  F(2kQLe  sin*  ^ ) 

Fl,/2(*oic(l  -  \/l  -  exp(-jifcosc)  „ 


rsin  p  cos  ( 


(3.19) 


3.4  Numerical  results  and  discussion. 


The  corner  diffracted  field  described  by  eqs.  (3.11a)-(3.11b)  does  not  satisfy 
either  the  boundary  conditions  on  the  surface  or  reciprocity  with  respect  to  the 
azimuthal  coordinate.  Both  of  these  inadequacies  are  inherent  with  the  Physical 
Optics  approximation  of  the  diffracted  wave.  Besides,  the  asymptotic  evalution 
of  the  radiation  integral  required  that  the  observation  point  was  not  close  to  the 
surface  of  the  semi-infinite  wedge,  where  the  direct  dipole  field  (edge  wave)  is 
dominant .  However,  the  Physical  Optics  solution  based  on  the  true  currents  that 
would  flow  over  the  surface  of  the  wedge  yields  a  satisfactory  approximation  of  1h« 
scattered  field,  provided  that  no  diffracted  rays  from  edges  (1)  and  (2)  reach  th** 
observation  point.  The  latter  corresponds  to  the  case  where  01,2  +  l^i,2  >  s’-  A 

a  ^ 

remarkable  fact  is  that  the  vertex  diflrari^d  HHd  is  dominantly  ^)-directed,  which  is 
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physically  justifiable  since  the  current  is  accumulated  in  the  vicinity  of  the  guiding 
edge,  flowing  parallel  to  it.  Note  also  that  as  a  — »  ir  the  vertex  diffracted  field 
vanishes  and  the  total  solution  reduces  to  the  edge  wave  term  as  it  should. 


Figure  12:  The  square  plate  used  in  the  numerical  application  of  the  Physical 

Optics  solution. 

The  radiation  of  a  dipole  in  the  close  vicinity  of  a  flat  square  plate  (see  Fig.  12) 
is  examined  in  Figs.  13-20.  The  dipole  is  located  at  a  distance  p' s  0.01  A  from 
the  center  of  one  of  the  edge  of  a  rectangular  plate  with  dimensions  4A  x  4A  as 
depicted  in  Fig.  12.  The  pattern  is  evaluated  in  the  elevation  plane  associated  with 
the  guiding  edge  fixed  reference  frame,  rentrred  at  the  point  of  projection  of  the 
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dipole  to  the  edge,  for  several  values  of  the  azimuthal  angle  ^  of  the  observation 
point.  In  particular,  the  radiated  field  is  approximated  by  the  direct  dipole  field 
(edge  wave)  multiplicatively  corrected  by  the  transition  dyadic  which  is  developed 

I 

ill  Chapter  IV,  and  the  superposition  of  two  corner  diffracted  rays  from  the  corners 
of  the  guiding  edge.  Other  possible  diffraction  mechanisms  which  may  contribute 
for  small  azimuthal  coordinates  of  the  observation  point  are  neglected  in  the  present 
analysis. 

The  predicted  radiation  patterns  are  compared  with  moment  method  results 
[9].  It  is  noted  that  moment  method  under  those  circumstances  should  not  be 
viewed  as  exact  since  its  surface  patch  monopoles  may  not  accurately  approximate 
the  singular  current  flowing  along  the  guiding  edge  as  well  as  the  strong  coupling 
between  the  dipole  and  that  edge.  As  an  independent  technique,  however,  it  yields 
comparable  results,  which  for  certain  aspects  of  observation  show  good  agreement 
with  the  Physical  Optics  approximation  developed  earlier. 

The  edge  wave  itself  and  the  edge  wave  corner  diffracted  field  clearly  dom¬ 
inate  the  total  0-diTtcitd  field  in  the  paraxial  region  of  the  guiding  edge.  Both 
mechanisms,  however,  decay  as  the  observer  is  removed  from  the  paraxial  region. 
The  noil-oscillatory  behavior  of  the  ^directed  field  amplitudes  (see,  for  instance. 
Figs.  19,  20)  also  reveals  that  the  corner  diffracted  contribution  to  that  compo¬ 
nent  of  the  field  appears  indeed  to  be  negligible.  Besides,  this  component  of  the 
radiated  wave  may  be  attributed  to  double  edge  diffractions  or  higher  order  terms 
in  the  asymptotic  expansion  of  the  radiation  integral  of  the  edge  wave  currents. 

The  comparability  of  the  moment  method  results  and  the  Physical  Optics  s** 
lution  deteriorates  at  small  azimuthal  angles.  However,  at  these  lower  cuts  other 
secondary  diffraction  mechanisms  (corner,  double  edxe  diffraction,  etc.)  may  con¬ 
tribute  significantly  to  the  radiation  pattern. 
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Pigvrv  13:  FIm  IWd  rcgioa  ^-dirvcted  w»ve  rMliftUd  by  Ibe  dipok  Pig«r*  14:  Fw  field  region  (l-direeted  wn»e  rndinted  by  tbe  dipole 

of  Pig.  12  nt  on  uimnthal  nagle  4  =  30*.  of  Pig.  12  ot  m  uimuihnl  Ugle  4  =  45*. 


Pig«re  17:  PWr  Md  Kgioa  ^-directed  w»»e  radikted  by  ib«  dipok  Figure  18:  Fur  Held  region  4' 
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CHAPTER  IV 


An  Equivalent  Current  Approach 


4.1  Introduction. 

It  was  shown  in  Chapter  III  that  the  surface  radiation  integral  can  be  asymp¬ 
totically  reduced  to  a  line  integral  along  the  terminating  edges  of  the  semi-infinite 
wedge.  In  other  words,  the  high  frequency  radiated  field  can  be  adequately  approx¬ 
imated  by  the  radiation  integral  of  equivalent  electric  and  magnetic  edge  currents 
flowing  along  these  terminating  edges.  In  the  Physical  Optics  approximation  de¬ 
veloped  earlier,  those  equivalent  edge  currents  asymptotically  incorporated  the 
physical  optics  induced  current  effect. 

Within  the  frame  of  the  Physical  Theory  of  Diffraction  [10,11,12]  the  scattered 
field  from  a  body  can  be  considered  as  the  superposition  of  the  waves  produced 
by  two  current  components  flowing  along  the  surface  of  the  body;  a  uniform  com¬ 
ponent  related  with  the  geometrical  optics  field  (physical  optics  currents)  and  a 
non-uniform  current  due  to  the  possible  curvature  or  discontinuities  in  the  curva¬ 
ture  of  the  surface  of  the  scatterer  (fringe  currents).  Hence,  in  the  particular  edge 
wave  vertex  diffraction  problem,  the  geometry  of  which  is  illustrated  in  Fig.  10,  the 
solution  can  be  improved  by  superimposing  to  the  uniform  currents  fringe  current' 
due  to  the  presence  of  the  terminating  edge  discontinuity.  In  the  context  of  our 
asymptotic  analysis  the  field  is  approximated  by  the  radiation  integral  of  the  sum 
of  two  types  of  edge  (both  electric  and  nisKnetic)  equivalent  currents:  the  first 
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set  corresponding  to  the  uniform  (physical  optics)  currents  and  the  second  incor¬ 
porating  the  fringe  current  effect.  The  asymptotic  reduction  of  the  edge  current 
radiation  integral  involves  a  stationary  phase  point  contribution  (edge  diffracted 

t 

wave)  together  with  an  end  point  effect.  The  former,  as  the  origin  of  Keller’s  cone 
of  diffracted  rays  moves  away  from  the  vertex,  and  in  the  region  that  is  illuminated 
by  the  edge  diffracted  field,  should  normally  reduce  to  the  UTD  solution,  whereas 
the  latter  is  interpreted  as  a  vertex  diffracted  wave. 

The  available  representations  of  the  edge  equivalent  currents  are  based  on 
the  canonical  solution  to  the  infinite  wedge  problem  or  its  asymptotic  reductions 
(Keller’s  Theory,  UTD),  which  clearly  fail  in  the  vicinity  of  the  the  vertex.  How¬ 
ever,  the  uniform  asymptotic  representation  of  the  total  field  obtained  from  the 
truncated  equivalent  currents  radiation  integral  presents  a  reasonable  qualitative 
behavior  and  shows  good  agreement  with  moment  method  and  measured  results. 

4.2  The  equivalent  edge  currents  concept. 

The  concept  of  the  equivalent  edge  currents  in  the  theory  of  diffraction,  based 
on  Sommerfeld’s  half  plane  solution,  was  initially  developed  by  Millar  [13,14,15] 
and  Clemmow  |]6|  as  a  simple  method  for  deriving  asymptotic  approximations 
of  the  electromagnetic  wave  diffracted  by  large  apertures  in  perfectly  conducting 
screens.  Later,  the  same  concept  was  formalised  by  Ryan  and  Peters  [17]  and 
employed  to  correct  the  singular  GTD  fields  in  the  caustic  regions.  A  UTD  repre¬ 
sentation  of  the  Ryan  and  Peters  equivalent  edge  currents  for  plane  wave  incidence 
on  a  perfectly  conducting  wedge  and  for  observation  points  sufficiently  far  from 
the  edge  of  the  wedge,  may  be  written  in  terms  of  the  diffraction  angle  i3q.  the 
distance  s  of  the  observation  points  from  the  point  Q  of  the  wedge  at  which  the 
currents  are  evaluated,  and  the  asimuthal  r«.»ordinales  of  the  diffracted  and 
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incident  ray  with  respect  to  the  edge  fixed  reference  frame,  as  follows: 


^HQ)  =  B‘(Q)  ■  (4.1.) 

A/'»(Q)  =  -^H'(Q)./m(.;*,#';A,)  (4.1b) 

with 


+  <-ol  F|to.«-{*  -  *')) 

t{coi(!!^^^)f(W(*  +  *')I 

+  “t  ~  *'*)  F|lto*.-(*  +  *'))}]  ,  (4.1c) 


0^(7)  =  2cos^ 


(4.1d) 


being  the  integers  that  most  nearly  satisfy  the  equations 


2rrn7V*  -  7  =  ±7r  (4.1e) 

and  E'(Q)i  ff*{Q)  denoting  the  incident  at  Q  electric  and  magnetic  field,  respec¬ 
tively.  As  before,  the  parameter  n  (wedge  number)  is  related  with  the  wedge  angle 
(WA)  via  ( VVA)  =  (2  —  n)7r.  Note,  also,  that  the  edge  transition  functions  F(-)  are 
not  normally  associated  with  the  Ryan  and  Peters  form.  It  should  be  emphasized 
that  Ryan  and  Peters  equivalent  currents  are  valid  for  directions  of  observation 
lying  on  Keller’s  cone  of  diffracted  rays.  Knott  and  Senior  [18]  suggested  a  gener¬ 
alization  of  the  expressions  of  the  equivalent  currents  in  order  to  include  arbitrary 
directions  of  incidence  and  observation  by  deforming  the  term  sin^  appearing  in 
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the  denominator  of  eqs.  (4.1a)-(4.1e)  to  the  product  sin/^sin/S'  with  ^  denot¬ 
ing  the  angle  of  observation  and  incidence  respectively  in  a  predefined  edge  fixed 
coordinate  system.  But  their  proposal  is  merely  a  postulate,  since  it  is  bued  on 
intuitive  symmentry  considerations  rather  than  any  mathematical  derivation. 

Recently,  Michaeli  [19,20,21,22]  proposed  explicit  expressions  for  the  equiva¬ 
lent  edge  currents  deduced  from  the  exact  solution  of  the  canonical  wedge  problem 
and  valid  for  arbitrary  directions  of  illumination  and  observation.  His  currents  are 
related  with  true  currents  flowing  along  the  wedge  surface,  since  their  derivation  is 
based  on  the  asymptotic  relationship  between  the  Physical  Theory  of  Diffraction 
surface  radiation  integrsJ  and  the  equivalent  current  line  integral. 

Equivalent  currents,  including  higher  order  interaction  effects,  were  utilized  by 
Sikta  [4]  for  the  computation  of  the  off  principal  plane  RCS  of  flat  plate  structures. 
Although  Sikta  used  a  modified  version  of  the  currents  of  eqs.  (4.1a)-(4.1e),  which 
is  conceptually  empirical  his  results  are  comparable  with  those  presented  later 
in  (22)  where  the  more  rigorous  Michaeli’s  edge  currents  were  used.  As  pointed 
out  in  Chapter  I,  Burnside  and  Pathak  |3]  employed  the  radiation  integral  of 
the  truncated  equivalent  currents  to  asymptotically  identify  a  corner  diffraction 
coefficient  after  an  empirical  modification  of  the  final  result.  A  similar  approach  is 
adopted  in  this  chapter,  using  two  types  of  equivalent  currents,  to  develop  an  edge 
wave  vertex  diffraction  coefficient  which  incorporates  the  singular  behavior  of  the 
edge  wave  in  the  vicinity  of  the  guiding  edge  and  an  edge  diffraction  coefficient 
valid  as  the  origin  of  Keller’s  cone  approaches  the  vertex  point. 
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4.2.1  Generalized  equivalent  currents.  The  half  plane  case. 
Non  grazing  incidence 


A  plane  wave  illuminates  a  half  plane  as  illustrated  in  Fig.  2.  The  problem  as 
already  pointed  out  is  essentially  two  dimensional.  The  total  field  can  be  expressed 
in  terms  of  the  vector  potentials 

y4  =  fV’"  ;  F  =  iV’^  (4.2) 

via  the  following  equations: 


E  = -V  y.  F  -  ju)QfiQ{l kQ^V\7)- A  (4.3a) 

H  =  V  X  A  -  jwQtoil  +  k^^W)  •  F  (4.3b) 


where  I  is  the  unit  dyadic.  Furthermore,  the  total  field  can  be  viewed  as  a  field 
created  by  fictitious  equivalent  line  sources  flowing  along  the  edge  of  the  half  plane, 
so  that  one  may  write 


r(i1=  r 

/-oo 

(4.4) 

v^(i*)=  r  M^{s;0\<t>'-,z')Go{p,z\z')dz’ 

i-oo 

(4.5) 

in  which 


exp{-jfco\/(z  -  z*)2  + 
Air^(z  ~  z')2  + 


(4.6) 


is  the  free  space  scalar  Green’s  function,  s  is  the  vector  pointing  to  the  observer 
in  an  arbitrarily  chosen  reference  frame,  0'  is  the  angle  between  the  incident  ray 
and  the  wedge  and  the  superscript  i  indicates  the  connection  with  the  total  field. 
It  is  also  reminded  that  the  p  and  ^  -directed  components  of  the  electric  and 
magnetic  field  can  be  expressed  in  terms  of  its  f-directed  components  as  indicated 
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mathematically  by  cqs.  (2.1a),  (2.1b).  Thus,  from  eqs.  (4.2)-(4.6),  one  obtains  for 
the  i-directed  components 

y — 00 

^2|/‘(i-i^',*';:')Co((>,z|V)|)i2’  (4.7.) 

N,(i)  =  -juixo 

y— oo 

+*=0*  (4.7b) 

which  can  be  employed  for  a  unique  solution  for  7^,  Af*.  Noting  the  convolution 
type  of  the  integrals  involved  in  the  integral  equations  (4.7a),  (4.7b)  the  latter 
can  be  formally  solved  with  the  aid  of  the  Fourier  transform  with  respect  to  the 
variable  z  (the  z  coordinate  of  the  observation  point).  It  is  also  remarked  that 
the  principal  value  of  the  integrals  involved  in  cqs.  (4.7a),  (4.7b)  with  respect  to 
the  arbitrarily  chosen  origin  on  the  l-axis  should  be  considered,  along  with  the 
definition  of  the  Fourier  transform  (23).  Without  presenting  the  details  of  the 
derivation  and  restricting  our  attention  to  the  diffracted  term  of  the  total  field, 
which  is  physically  related  with  currents  flowing  in  the  vicinity  of  the  edge,  one 
obtains  the  following  expressions  for  the  edge  currents: 

I  _  2exp(-yir/4) 

_ expf-yibossinjgsin/?') _ 

y $  sin  0  tin  0'  kos  sin  0  sin  0' ) 

nF*(z')^zifis;0\^’;z') 
ZoH‘(z').im/(s-^',<^'-,z') 
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Now,  the  superscript  /  indicates  relation  with  the  diffracted  field  and  the  distance 
a  is  related  with  the  radial  coordinate  p  via  p  —  a  Moreover, 

F(2fco<  *in  0  sin  cos^ 
cos 

Note  that  the  same  result  can  be  obtained  by  merely  dividing  the  diffracted  field 
with  tlie  two  dimensional  free  space  Green’s  function.  For  large  arguments  of  the 
Hankel  function  and  if  the  origin  of  the  i*axis  coincides  with  the  diffraction  point 
Q*  the  equivalent  edge  currents  reduce  to  Ryan  and  Peters  expressions  presented  by 
eqs.  (4.1a)-(4.1e).  One  can  modify  the  above  expressions  for  arbitrary  directions  of 
incidence  and  observation  and  spherical  wave  incidence  by  substituting  the  factor 
sin^/?'  appearing  in  the  denominator  of  eq.  (4.8b)  with  sin  ^  sin  and  a  with  tne 
distance  parameter  L  —  aa* /(a-^-  s'),  where  s'  denoting  the  distance  of  the  source 
from  the  point  Q  at  which  the  equivalent  edge  currents  are  evaluated.  Then,  for 
large  values  of  the  parameter  !>,  one  obtains 


where  now 


sin  0  sin  jV 


F\2koL  sin  0  sin  0'  cos*  ^ ^  ] 
cos  (^) 


F\2k^L  sin  0  sin  0'  cos*  ^  j  ] 


(4.9b) 


and  i  is  again  the  unit  vector  tangent  to  tb»'  edge  nt  Q. 
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The  perturbation  of  the  original  expressions  (4.8a),  (4.8b)  is  based  on  the 
same  arguments  stated  by  Knott  and  Senior  and  has  no  mathematical  justifica¬ 
tion.  Evidently,'  the  above  representation  can  be  generalized  for  the  wedge  rase. 
The  final  result  is  merely  a  symmetric  perturbation  of  Ryan  and  Peters  edge  cur¬ 
rents.  It  should  be  emphasized  that  asymptotic  evaluation  of  the  radiation  in- 
-  tegral  of  the  edge  currents  described  by  eqs.  (4.9a)-(4.9b)  over  a  finite  (possibly 
curved)  edge  yields  always  a  reciprocal  result  and  the  correct  field  for  observation 
directions  in  the  Keller’s  cone  of  diffracted  rays.  Besides,  the  more  general  expres¬ 
sions  (4.8a),  (4.8b)  can  be  used  for  corrections  in  caustic  regions  in  the  near  field 
of  the  wedge  where  the  GTD  related  Ryan  and  Peters  equivalent  currents  fail. 

Grazing  incidence 

An  interesting  case  arises  when  the  incident  plane  wave  grazes  the  edge  of 
the  half  plane  (/?'  -*  0).  The  equivalent  edge  currents  that  would  support  the 
edge  guided  wave  given  by  eqs.  (2.13),  (2.14)  can  be  derived  following  a  similar 
procedure.  However,  in  this  case,  the  integral  equations  (4.7a),  (4.7b)  as  — *  0 
yield  an  indeterminate  form  which  implies  an  infinity  of  solutions  for  the  unknown 
edge  currents.  This  ambiguity  originates  from  the  fact  that  grazing  plane  wave 
incidence  is  essentially  a  theoretical  idealization.  Nevertheless,  based  on  the  dipole 
excitation  problem  examined  in  Section  2.3,  one  can  obtain  expressions  for  the 
equivalent  currents  which  are  physically  reasonable.  Beginning  with  the  integral 
equations: 

(4.10a) 
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(4.10b) 


as  i3'  — »  0,  ontf  may  postulate 


r^’(s;0\<l>’;z) 


8exp(j>/4)  y/pSo{^\ii>') 
H^^^kop  Bin  0>) 


txp{-jkoz) 


lo8in(^/2) 
-  C05(^/2) 


(4.11) 


where 


c  /  y  ^00'  «"(^V2)  +  Eq^,  co8(^72) 

)  =  - 7== - 


(4.12) 


is  an  excitation  related  factor. 

4.3  Michaeli’s  edge  currents  in  an  oblique  edge  fixed  frame. 

For  the  general  wedge  case  one  obtains  the  following  equations  for  the  total 
edge  magnetic  and  electric  equivalent  currents  excited  at  the  point  Q  of  the  edge 
of  the  wedge  [19]: 


HQ)  =  -E'{Q)i 


JYq  sin(^7n)  f 
nto  sin^  0'  \ 


cos|(x  -  7i)/n]  -  coB{<f>'/n) 


— 1  -  SHQ)  •  <  — — — 

y/n)j  nko  Bin0' 


co8[(jr  -  72)/”]  +  cos(^7*») J  nkot 

{sin[(ir  -  7i)/n]  pj  cot  ^  —  cot  0  cos  ^ 

sin  71  cos|(x  -  7i)/n]  -  co8(^7n) 
8in[(ir  -  72)/n]  P2  cot  0'  —  cot^cos(nx  —  ^)  1 
sin  72  cos|(x  -  72)/n]  +  cob(4>' /n)  ] 


(4.13a) 


M(Q)  =  H\Q)i 


2j^o 

nibo  sin  0  sin  0' 
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3 


im((w  -7i)/n]  sin^ 


tin  12  {fOt[(ir  -  12)/»«1  +  cot(^7n)} 
where  the  complex  angles  71^2  defined  by  the  equation 


Figure  21:  Oblique  edge  fixed  coordinate  system  for  the  definition  of  Michaeli't 

equivalent  currents. 


71.2  =  -iln(Mi,j  +  Jmi,2  -  1) 


(4.14 


along  with  the  following  choice  of  the  branch  of  the  square  root  in  eq.  (4.14): 


|y//i*-  1|  ;  >  1 

;Vp2-l|  ;  -!</,  <1 


(4.1.5) 
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The  parameter  fj  depends  on  the  coordinates  (3,  0  and  <i/  of  the  observer 

and  the  source,  respectively,  in  terms  of  the  edge  fixed  reference  frame  centered 

at  Q  as  well  as  the  choice  of  the  edge  fixed  system  of  coordinates  itself,  i.e.,  the 

angle  q  between  the  unit  vector  t  tangent  to  the  edge  at  Q  and  the  unit  vector 

d  tangent  to  face  1  of  the  wedge  at  Q,  as  shown  in  Fig.  21.  The  latter  implies 

a  nonuniqueness  in  the  presentation  of  the  Michaeli’s  equivalent  currents.  For  a 

general  a  angle,  is  found  equal  to 

sin  Q  sin  l3  cos  4>  +  cos  a  (cos  —  cos  /3') 
sin  a  sin  /?' 

For  the  definition  of  ^<2  which  is  associated  with  face  2,  should  be  substituted  by 
niT  -  <f>  and  /?,  /?'  by  tt  -  /5,  rr  -  respectively.  If  the  classical  edge  fixed  system 
of  coordinates  is  chosen,  then  a  =  ir/2  and  /t]  becomes 

/il(ir/2)  =  sin  ^3  cos  <^/ sin (4.17) 


(4.16) 


On  the  other  hand,  if  the  d-axis  is  chosen  parallel  to  the  diffracted  ray  that 
grazes  face  1  surface,  as  suggested  by  Michael!  in  [20],  he.,  when  a  ~  from 
eq.  (4.15)  one  obtains 


Ml(^') 


sin  /?  sin  cos  <f>  -f  cos  fS'  (cos  0  —  cos  0' ) 


(4.18) 


For  a  more  systematic  approach  to  the  approximation  of  the  edge  wave  vertex 
diffraction  mechanism,  it  is  convenient  to  extract  the  physical  optics  components 
from  the  total  equivalent  edge  currents  presented  by  Michaeli.  In  accordance  with 
[24],  the  total  currents  given  by  eqs.  (4.13a),  {4.13b)  can  be  split  into  the  Physical 
Optics  current  components  A/P®,  /P®  and  the  fringe  current  components 


with 


/(<?)  =  I’”[Q)  +  I>(Q)  ;  M(Q)  =  +  M>(Q)  (4.1P) 
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Specifically,  the  Physical  Optics  components  associated  with  face  1  of  the  wedge 
are  given  by  the  following  equations: 


/r«?)  =• 


2j  U{n  - 


eHq)  •  / 


sin^' 
Zo  sin  0' 


kosml3*{ni  +  cos^') 

-i}*(Q)  •  i {coi 0^ cos +  C0t/?C08^)| 

2jZo  8in^f^(ir  -  ^') 


(4.20a) 

(4.20b) 


ka  sin/?  6in/?'(pi  +  cos^') 
while  the  fringe  currents  can  be  derived  from  eqs.  (4.13a),  (4.13b),  (4.19) 
and  (4.20a),  (4.20b). 

4.4  Vertex  dififraction  of  an  edge  wave  excited  by  a  dipole. 


An  excitation  dipole  is  placed  in  the  close  vicinity  of  the  edge  of  a  semi-infinite 
wedge  and  at  a  distance  (with  kos'^  >>  1)  from  its  vertex.  The  interaction 
between  the  dipole  and  tl’?  edge  produces  an  edge  guided  wave  which,  sufficiently 
far  from  the  dipole,  can  be  treated  as  ray  optical.  This  ray  optical  ^dge  wave 
is  diffracted  in  accordance  with  the  generalized  Fermat’s  principle.  In  our  case, 
despite  the  direct  field  which  is  essentially  the  edge  wave  suitably  modified  for  a 
semi-infinite  structure,  the  observer  is  reached  by  a  vertex  and,  possibly,  two  edge 
diffracted  rays  emanating  from  the  vertex  point  and  the  terminating  edges 
points  Qj,  Q2,  respectively,  where  Q|,  Q2  are  the  origins  of  the  classical  Keller’s 
cone  of  diffracted  rays.  The  observer  is  in  the  illuminated  or  shadowed  region  of 
the  edge  wave  edge  diffracted  fields  from  the  points  Qp  Q2’  ^2  points 

of  the  semi-infinite  edges  (1),  (2)  or  points  of  their  extensions,  respectively. 

A  uniform  solution  is  investigated  herein  using  the  equivalent  currents  de¬ 
scribed  earlier  in  this  chapter  and  exvited  by  the  edge  wave  impinging  on  the 
terminating  edges  (1)  and  (2).  The  “uniformity”  of  the  solution  requires  the  con¬ 
tinuity  of  the  total  field  at  the  shad*  v:  boundary  of  the  direct  field  (edge  wave) 
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and  the  boundaries  of  the  regions  illuminated  by  the  edge  diffracted  rays  from 
C?2*  Furthermore,  the  solution  should  remain  finite  at  the  extension  of  the 
guiding  edge.  One  may  write  for  the  total  field 


E*  =  •  r  +  Ef  )  +  E\  U{tQt )  +  E^ 


(4.21) 


where  is  the  field  produced  by  the  dipole  in  the  presence  of  an  infinite  wedge, 
7*  is  a  heuristic  transition  dyadic  which  assures  the  continuity  across  its  shadow 
boundaries  and  its  finiteness  along  the  extension  of  the  guiding  edge,  Ef  2  are  the 
edge  diffracted  fields  with  ^  denoting  the  i  coordinate  of  the  points  Q\  2  in 
terms  of  the  edges  (1)  and  (2)  fixed  coordinate  systems  depicted  in  Fig.  10,  and 
E^  is  the  vertex  diffracted  wave.  It  is  reminded  that  V{‘)  represents  the  unit  step 
function  (  V(x)  =  1  if  *  >  0,  U{x)  =  0  if  x  <  0). 

As  before,  our  attention  is  restricted  to  face  1  of  the  wedge.  The  equivalent 
currents  presented  earlier  for  grazing  edge  wave  incidence  at  the  point  Q\  of  the 
edge  ( 1 )  may  be  rewritten  as 


/l(Ql)  = 


2j 

fco  2 


(4.22a) 


*0 - 2 - 


(4.22b) 


where  E''‘'(Q])  is  the  field  produced  by  the  dipole  radiating  in  the  presence  of 
an  infinite  wedge  evaluated  at  the  point  Qi,  fj  is  the  unit  vector  tangent  to  the 
edge  at  Qi  and  tj,  m]  are  known  and,  in  genera],  slowly  varying  functions  of  the 
observer's  location  in  the  vertex  fixed  coordinate  system  and  the  distance  of  Qi 
from  the  tip  of  the  trihedron.  The  wedge  number  U]  is  related  now  with  the 
wedge  formed  by  face  1  and  the  plane  surface  defined  by  the  edges  (1)  and  (2). 
Note  also  that  a  factor  of  1/2  has  been  introduced  due  to  grazing  incidence.  It  can 
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be  easily  shown  using  eq.  (2.19)  and  the  relationship 


which  is  valid  for  ray  optical  waves,  that  ‘ 

^  "(Qi  )i>  =  -y*i(  «■>■'  a,  CW)  aUp',*')  ,y-i 

**1 

with 

ri  =  ysj*  +  <1  -  2fi  Sj  rosai  . 


(4.23) 

(4.24) 


As  before, 

K{p\  +  Pji  fo*  t'^) 

and  the  constant  C{u)  is  defined  by  eq.  (2.17). 

It  is  presumed  that  the  edge  and  vertex  diffracted  edge  wave  associated  with 
face  1  of  the  wedge  can  be  approximated  by  the  radiation  integral  of  the  equivalent 
currents  /i(^i),  Mi{Q\)  flowing  along  the  edge  (1),  which  in  the  Fresnel  or  the 
Fraunhofer  region  of  the  edge  can  be  explicitly  written  as 


■  j  (A  X  ii  X  + A  X  fimi(rc;(i)) 


.*/-!  exp{-i*o(^  +  n)) 


(4.25) 


with  R  =  llsc  ~  fifill  —  <1^1  )/^ 

For  large  values  of  the  parameter  fco  the  asymptotic  evalution  of  the  radiation 
integral  reveals  a  vertex  contribution  and,  also,  an  edge  diffracted  term,  if 


COSO]  (1  -  cos/J)  +  sinoi  sin/?  cos^  >  0  . 


i.e.,  if  Q\  lies  on  the  edge  (1)  itself  rather  than  on  its  extension.  These  two  term!' 
are  examined  separately  in  Sections  4.4.1  and  4.4.2  respectively. 


4.4.1  Edge  wave  vertex  diflfracted  field. 

Interpreting  the  end  point  contribution  to  the  radiation  integral  as  a  vertex 
diffraction  term, 'one  obtains  the  approximate  relationship 

s  jkS Zt  .in'' at  +  *cl> 

4ir»c  •c' 

•^(rc)/£.(lto)  (4.26) 

where  we  have  substituted 


^l(«c)  =  ic  X  ic  X  f'l  ti(/c;0)  +  ic  X  <j  mi(sc;0) 

and  denotes  the  end  point  contribution  to  the  integral 

^  exp{-ifco(^  +  *‘o)} • 

Referring  to  eq.  (3.8)  one  writes 

r(i/)tj‘'exp(>i/ir/2)«p(->M»c  +  *i)) 

_ 

[sina]  sin/Jcos^  +  co8ai(l  —  cos^)]^ 


(4.27) 


(4.28) 


in  which  the  edge  wave  transition  function  has  been  defined  by  eq.  (3.13).  The 
vertex  diffracted  field  associated  with  face  1  of  the  wedge  becomes 

£51.-,)  «  «pO(..  +  i)»/2)r(r)C(r)Z,4'(p',/)Stf^f-i^ 

4ksc*c 


irii/(fcooi(<c;<c)] 

(cot  ai(l  —  cos/?)  +  sin/?  cos 


(4.29) 


Michaeli’s  equivalent  currents 


It  has  been  pointed  out  that  Michaeli's  equivalent  line  sources  can  be  derived 
from  the  asymptotic  reduction  of  the  surface  radiation  integral  of  the  true  surface* 
currents  evaluated  from  the  canoiii«-»l  s«*ln«ion  miuI  flowing  over  the  faces  of  the 
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wedge  and  can  be  split  into  geometrical  optics  and  fringe  current  components. 
Moreover,  when  these  currents  are  truncated  to  reveal  an  end  point  effect,  one 
should  be  cautious  about  the  choice  of  the  parameter  fi  involved  in  the  expressions 
of  these  currents,  i.e.,  the  choice  of  the  edge  fixed  coordinate  system,  so  that  the 
terminated  equivalent  edge  currents  represent  correctly  the  end  point  effect  of  the 
truncation  by  the  edges  of  the  trihedron  true  surface  induced  currents.  In  this 
context,  a  correct  choice  of  the  parameter  fii  for  the  Physical  Optics  component 
of  the  equivalent  currents,  using  eq.  (4.16),  is  the  following: 


/^r  =  = 


sinaj  sin  01  cos^j  +  cosoi(cos/3i  +  cosoi) 


sin^  oi 


(4.30) 


so  that  the  unit  vector  &  in  Fig.  21  is  parallel  to  the  guiding  edge,  whereas  for  the 
fringe  currents  the  proper  choice  is 


f  ,  .  sinoi  sin/?i  cos^l  —  co5ai(co8/9i  +  cosai ) 

(il  =  Hl{ir  -  ai)  = - T-s -  (4.31) 

sin^oi 

and,  now,  &  is  parallel  to  the  edge  diffracted  ray  from  that  grazes  the  plane  of 
face  1.  In  the  above  equations  <f>i  are  the  elevation  and  azimuthal  angle  of  the 
observation  point  in  terms  of  the  'dge  (1)  fixed  coordinate  system  centered  at  Q^. 
The  previous  discussion  also  suggests  the  decomposition 


G\{sc)  =  G\^{sc)  +  G]'f{»c) 


(4.32) 


where  the  superscripts  po  and  /  imply  connection  with  the  physical  optics  and 
fringe  current  components,  respectively.  Not  surprisingly,  using  eq.  (4.30)  for  the 
definition  of  the  parameter  /»  in  the  eqs.  ('1.2nn)-(^.20b)  of  the  physirsd  optirs 
equivalent  edge  currents  and  after  some  tedious  manipulation,  which,  however, 
involves  only  elementary  operations,  it  can  be  shown  that: 

=  -0  cot  I  (4.33) 
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80  ih&t  the  vertex  diffracted  field  assumes  the  approximation 


(4.34) 


in  which  £^'^{jc)  >s  the  field  given  by  eq.  (3.10)  or  eqs.  (3.11a),  (3.11b),  while  the 
fringe  currents  related  wave  is  equal  to 

4irscSc‘'  2' 

.  ^i'(‘c)njkcci(rcy,)] 

(cot  Qi(l  —  COS  0)  +  sin  0  cos 

(cot  02(1  —  cos/?)  +  8in^cos(njr  —  ^))*' 

The  vector  functions  ^'^(sc)  (I  =  1,2)  are  related  with  the  Michaeli's  equivalent 


currents  at  via  the  equations: 


Jc)  =  i{ {Qc)ic  X  ie  X  +  m{ {Qc)ie  x  ii 


(4.36) 


where 


*/«?'■) 


um 

»”i(Q‘') 


(4.37) 


and  ij(Q^),  »n(((?^),  *|”(Q‘^)i  »n|”(Q^)  are  described  explicitly  by  the  expressions 


iMU - 

njsino/  I  siu 


cot/?/ cos 8in((y  -7/i)/n/] 
711  1 -ros((ir -7/i)/n/] 


^ify  cot  c^  +  cot  01  co8(n/ir  -  ^/)  sin((7r  -  7/2  )/n/]  | 

•iin/2  1  -  ‘)/2)/”l]  j 


(4.38) 


50 


and 


1  /  linKtr  -  7il)/w/l 

n|  sin  aj  »h\0i  \ 811171,  1  “  «o»l(»  “  7ll)/"ll 
8iii(n|ir  -  ^i)  8in((ir  -  712)/"/)  \ 

tin 712  1  +  cos{(ir  -  7/2)/*»l]  / 


(4.39) 


iPOinc\  _  “  cot/9|  co8^ 


*n^') = 

= “T 


tina,(l  +/if) 
8in^/ 


(4.40) 


(4.41) 


sma/Bin/9/(l  +  M| ) 

Tli«  index  /  assumes  the  values  1  and  2  for  the  edges  (1)  and  (2)  of  the  trihedron, 
respectively.  The  elevation  and  azimuthal  angles  ^,,2*  ^1,3  observer  in  terms 
of  the  coordinate  systems  associated  with  the  edges  (1)  and  (2)  and  centered  at 
the  tip  are  related  with  the  principal  spherical  coordinates  /9,  ^  via  the  easily 
verified  equations 


cos^l  =  sinai8in/9co8^-cosaico8/9 ; 
9\n0i  —  1^1  -  C08^/?l| 


(4.42) 


cos  02  —  ®2  **“  cos(nir  —  ^)  —  cos  02  cos  0  ; 

tin  ^2  =  Iv^l  -  cos^^il 


(4.43) 


COS^l  = 


■in^l  = 


nndriii^ 

SIB01 


(4.44) 


,  eosoi sin ^co»(ii»-4Ksino2 «•©•/?  . 

cos^2= - ^ - riadj 

In  nccordnnre  with  eq.  (4.31)  (h«  pnriuiiftrr.  p/j.  p/,  *ltt>uld  be  rhoeen  equkl  to 
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J  - 

sin  o/  Bin  01  cos  —  cos  q/(cos^I  +  cos  0/) 

(4.46) 

Nl  - 

sin'  0/ 

J  - 

sin  ai  sin  0i  cos(n/jr  —  ^i)  —  cos  a/(cos  0i  -f  cos  a/) 

(4.47) 

h2  - 

sin'  0| 

t 

The  above  expressions  simplify  considerably  for  the  case  of  the  plane  angular  sector, 
with  aj  =02  =  o  ,ai  =  02  =  «■  Without  presenting  the  details  of  the  derivation 
which  again  involves  only  elementary  manipulations,  for  the  field  related  with  the 
fringe  currents  one  derives 


a  a)  i^i/2l*0«*(*e;  4)1 

F{2koLe  sin*  ^ ~  ~  ^  ~ " 

2  ^cot  q(1  —  cos  fi)  +  sin  fi  cos  ^ 


where  the  function  G^'f(»c)  i«  given  by 


1 


cos(2o)  +  cobS 


(4.48a) 


(4.48b) 


with 


cos^  =  sin(2Q)sin/3cos^  -  cos(2a)cos^ 


(4.48c) 


*0(0-,^)  =  ~( cos  a  +  cos /?i)  (sin  a  cos /3  cos  ^  + COSO  sin /9) 

•(sinocos^l(l  -  cosocos^l)  -  eos*a8in/?iI 
■f  sin*osin^sin^l 


(4.48d) 


=  sinosin^(cosa  +  co8/.?i) 

•(sin  o  cos  1  -  cos  o  cos  /^j )  -  cos*  a  sin 0i\ 
+  sin^  o  sin  (sin  o  cos  /.t  cos  +  cos  o  sin /}) 


(4.48e) 
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The  vertex  diffracted  field  derived  with  Michaeli’s  equivalent  currents  has  a 
clear  physical  interpretation  in  that  it  is  essentially  the  field  produced  by  the  super¬ 
position  of  the  truncated  physical  optics  currents  that  would  flow  over  an  infinite 
wedge  and  the  fringe  currents  excited  by  the  terminated  edges  for  edge  wave  grazing 
incidence.  However,  it  does  not  satisfy  the  boundary  conditions  on  the  trihedron 
surface.  Although  the  latter  may  be  a  minor  concern  for  Fraunhofer  region  obser¬ 
vations,  as  far  as  the  earlier  derived  equations  serve  as  a  good  approximation  of  the 
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region  observations,  as  far  as  the  earlier  derived  equations  serve  as  a  good  approxi¬ 
mation  of  the  vertex  diffracted  edge  wave,  it  would  be  interesting  to  obtain  a  vertex 
diffraction  coefficient  which  would  satisfy  the  appropriate  boundary  conditions  on 
the  surface  of  the  trihedron.  For  this  purpose  one  may  introduce  the  generalised 
equivalent  currents  described  for  the  half  plane  case  by  eqs.  (4.9a)-(4.9b).  The 
approach  is  essentially  empirical  since  the  generalised  equivalent  currents  are  not 
necessarily  related  with  true  surface  currents,  but  merely  is  a  symmetric  perturba¬ 
tion  of  the  edge  currents  that  would  produce  the  correct  edge  diffracted  field  irouli 
an  infinite  edge. 

For  simplicity  the  plane  right  angular  sector  rather  than  the  complex  trihedron 
problem  is  examined.  The  total  vertex  diffracted  field  is  related  with  the  end 
point  contribution  to  the  radiation  integral  of  two  equivalent  line  sources:  the  first 
associated  with  the  guiding  edge  (edge  (0))  which  is  grased  by  the  dipole  excited 
field  and  can  be  derived  by  a  heuristic  modification  of  the  currents  of  eq.  (4.11), 
and  the  second  associated  with  the  terminated  edge  (edge  (1))  and  evaluated  from 
eqs.  (4.9a)-(4.9b)  for  edge  wave  incidence.  The  edge  wave  equivalent  currents  given 
by  eq.  (4.11 )  are  clearly  valid  only  for  plane  wave  incidence  which  is  also  evident 
in  their  dependence  rather  than  the  dependence  of  Ryan  and  Peters  as 
well  as  Michaeli  equivalent  currents.  The  empirical  modification  concerned  with 
the  edge  wave  dipole  excitation  problem  involves  a  multiplicative  introduction  of 
a  "reflection  coefficient”  similar  to  that  suggested  by  Sikta  in  (4],  namely 

4  yfirkQs'c 

as  well  as  the  multiplication  by  the  edge  wave  transition  function 


(4.521 


which  assures  the  finiteness  of  the  total  solution  as  0  —*  0.  The  corner  diffracted 
field  associated  with  the  guiding  edge  is  then  equal  to 

j  sinj  +  E’^jCT)  ,mi 


E^{3c)  W 


y/2irkQ 


yjsin0  sin /S' 


.  1  <f>^F{2koLcsin^ 
ip  ««  cos  -  ) - ^ ^ 


sin ' 


(4.53) 


Note  that  the  field  exhibits  a  similar  qualitative  behavior  with  the  edge  wave  corner 
diffracted  field  derived  in  [4]  from  a  limiting  manipulation  of  the  corner  diffraction 


coefficient.  The  term  associated  with  the  terminated  edge  is  written 

1  )  .in  ^  +  4((y)  cos  i 


\/5  irfco 

•(/3  sin  ^  ^  cos  /3  cos 


F(2fco£c  sin^  fii  cos^  ^ ) 
sin  01  cos  ^ 


F^^^2[^QLe{l  -  \/l  -  sin^  0 cos^  <^)) 


•v/sin  0  cos  (j> 


(4.54) 


The  total  vertex  diffracted  field,  however,  cannot  be  simply  obtained  from  the 
superposition  of  Eq  and  E^,  because  the  effect  of  the  truncated  currents  may  be 
double  counted.  Therefore,  a  modification  of  the  final  expressions  is  necessary  in 
order  to  obtain  results  compatible  with  those  presented  earlier. 

Although  the  new  expressions  for  the  field  satisfy  the  appropriate  boundary 
conditions  on  the  half  plane,  they  exhibit  a  discontinuity  on  the  plane  of  the  sector 
as  the  boundary  ^  =  0  is  crossed,  i.e.,  the  field  has  a  definite  finite  value  as 
0  0  on  the  =  it  half  plane,  while  it  vanishes  on  the  <j>  =  0  semi-plane.  This 

inadequacy  is  also  present  in  Sikta’s  expressions  of  the  edge  wave  corner  diffracted 
field.  Note  also  that  the  field  expressions  exhibit  a  ^directed  component  which 
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does  not  vanish,  as  it  should,  on  the  plane  of  the  plate.  The  same  problem  occurs 
in  the  expression  of  the  corner  diffraction  coefficient  and  is  coherent  with  the 
nature  of  the  truncated  equivalent  currents  which  are  derived  from  the  half  plane 

t 

solution.  It  can  be  possibly  overcome  by  incorporating  higher  order  interaction 
effects  between  the  currents  flowing  along  edge  (1)  and  the  adjacent  edge  (edge 
(0)),  which,  however,  appears  cumbersome. 


I -  2X  - 1 


Figure  22:  Angular  sector  geometry  for  the  comparison  of  the  generalized  and 

Michaeli's  equivalent  edge  currents. 

The  generalized  equivalent  current  field  associated  with  the  edge  (1)  (GEC 
(1))  is  compared  with  the  Michaeli’s  current  (MEC)  and  Physical  Optics  (PO) 
solution  in  Figs.  23  through  28  for  an  angular  sector  with  a  —  120^,  the  geometry 
of  which  is  shown  in  Fig.  22.  The  direct  dipole  field  is  not  incorporated  in  any 
of  the  above  calculations.  The  singularity  of  the  field  predicted  by  the  Physical 
Optics  solution  is  evident  in  Figs.  23  and  2^.  The  three  approaches  show  good 
agreement  for  small  values  of  the  elevati«»n  ans'l''  and  sufficiently  far  from  the 
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plane  of  the  angular  sector.  At  ^  =  180®  the  generalized  equivalent  current  field 
presents  a  non- vanishing  <^-directed  component,  which  is  not  physically  acceptable. 
However,  it  satisfies  the  boundary  conditions  on  the  sector  surface.  The  solutions 
diverge  at  the  vicinity  of  the  guiding  edge  (i.e.,  when  ^  — »  180®).  But  in  that 
region  the  total  field  is  dominated  by  the  singular  direct  dipole  field  (edge  wave). 
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Figwr*  25:  Fbr  Md  ^-dirrclrd  component  nt  mi  election  Migle  Figure  26:  Fur  field  ^directed  component  at  mi  eiemtion  Migle 

^  s  90*  for  the  angulM  eector  geometry  of  Fig.  22.  =  90®  for  the  MigulM  •eclor  geometry  of  Fig.  22. 
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4.4.2  Edge  wave  edge  diffracted  field. 

In  the  genera]  wedge  problem  the  point  Qj  lies  on  the  semi-infinite  edge  (1)  if 
rota|(l  —  cos/?)  +  sin/9ro8^  >  0  . 


The  angle  of  diffraction  /?j  in  this  circumstance  is  defined  by  the  equation 


(4.55) 


where 


I 


2  tq*  Sf  COS  ai 


(4.56) 


and  the  distance  tpe  of  the  point  Q|  from  the  tip  is  given  by 

^  sc4  «n(ai +/3i)  _ 

=  -f—- - ; - r— r  (4.5<) 

Sf.  sin  aj  +  Sc  wn  p\ 

The  non-uniform  representation  of  the  edge  diffracted  edge  wave  from  Q\ 
(non-uniform  with  respect  to  the  singular  behavior  of  the  impinging  edge  wave  as 
well  as  the  discontinuity  at  the  boundary  of  the  illuminated  by  the  diffracted  field 
region)  can  be  obtained  either  using  UTD  or,  equivalently,  evaluating  the  isolated 
stationary  phase  point  contribution  to  the  radiation  integral  of  the  edge  currents. 
(Note  that  both  types  of  equivalent  currents  utilized  earlier  for  the  approximation 
of  the  vertex  diffracted  field  result  in  the  Ryan  and  Peters  equivalent  currents  when 
evaluated  at  the  point  Qj ).  For  instance,  employing  UTD  the  edge  diffracted  edge 
wave  reads 


exp(~jx/4) 
2n  j  \/2irIo  sin 


+  Flk^Lia- 
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«l(«J  +  *1) 


exp(-ifcosi) 


(4.58a) 


where 


f  sin^/Sf 
*1  +  -1  * 


(4.58b) 


»\  =  yjiqt^  +  -  2<^«JeC08/3i  (4.58c) 

The  incident  edge  wave  can  be  derived  from  eq.  (2.19)  and  is  given  by 

fur  ccf )  =  J*o  ^0  C{U)  Al(p\  t')  «■!-'  oi  (4.59) 

Sj 

Note,  also,  that  a  factor  of  1/2  has  been  introduced  into  eq.  (4.58a)  due  to  grazing 
incidence  |25].  Clearly,  as  ^\{*{)  becomes  infinite.  The  singularity  of 

the  field  in  accordance  with  the  development  in  Appendix  A  can  be  compensated 
with  (he  multiplicative  introduction  of  (he  transition  function 


where 


/;.,(■)  =  (4.60) 

witli  th«  itar  4enoltiig  complex  conjogate.  The  large  parameter  Xe,  equals 

i.,  =  (4.61) 

Besides,  the  same  modification  factor  insures  the  continuity  of  the  total  diffrscled 
field  (edge  and  vertex  difracted  waves)  at  the  boundary  plane 


cos  o  1  ( 1  -  cos  /3 )  +  sin  o  |  sin  /3  ms  I) 


(4.62) 
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To  illustrate  the  resulting  continuity  of  the  diffracted  field  at  this  boundary  the 
particular  case  of  the  plane  angular  sector  {u  =  1/2;  ai  =  a)  is  examined.  In  this 
circumstance  tlie  proposed  uniform  approximation  of  the  edge  diffracted  field  may 
be  written  in  the  following  manner: 


^  =2 - r-. 


sin  Q  +  Se  sin 
in  a  ^  Scs'.si(si  +  Sj) 


exp{-ifco(*l  +  «i)} 


F{2koLesin^0\  cos^^)  Fi^koL^  cos"^  (^^))  . 


sin  01  cos  ^ 


y8in(o  +  /?i) 


(4.63) 


On  the  other  hand,  using  the  fact  that  Michaeli's  equivalent  currents  reduce  to 
those  of  Ryan  and  Peters  at  the  origin  of  Keller’s  cone  of  diffracted  rays,  or  alter¬ 
natively  employing  eqs.  (3.18),  (4.34),  (4.48a )-(4.48e)  (for  a  =  oii),  it  can  be  easily 


shown  that  the  corner  diffracted  field,  as  equals 

_  jZQA^p\<f>')  cxp{--jfco(^  4)}  F{2koLenn^ 


iir^vsma 


■  lim  - '—j===== - 

/?!  -•v-o  v'<‘os  o  +  cos  01 


(4.64) 


Then  based  on  the  definition  (3.13)  of  the  edge  wave  transition  function  and  the 
relationship  (4.63)  it  is  trivial  to  show  that  the  total  diffracted  field  remains  con¬ 
tinuous  as  the  boundary  of  the  region  illuminated  by  the  edge  diffracted  field  is 
crossed.  The  transition  function  i^_](*)  for  the  general  wedge  case  appears  as  a 
type  of  caustic  correction  factor  in  the  sense  that  it  compensates  the  singularity  of 
the  edge  diffracted  field  at  the  extension  of  the  guiding  edge.  However,  in  our  case 
the  singularity  arises  from  the  behavior  of  the  incident  field  on  the  edge  rather 
than  the  focusing  of  the  diffracted  rays  into  a  caustic. 
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4.4.3  A  transition  function  for  the  guiding  edge  diffracted  field. 

As  pointed  out  in  Section  2.5  the  rigorous  solution  of  the  radiation  of  a  dipole 
in  the  presence'  of  an  infinite  wedge  predicts  a  singular  field  at  the  edge  of  the 
wedge,  in  consistency  with  Meixner’s  edge  condition.  However,  this  singularity  of 
the  edge  wave  is  not  physically  acceptable  in  the  extension  of  a  semi-infinite  edge. 
For  overcoming  this  discrepancy  ihe  multiplicative  correction  of  the  edge  wave 
associated  with  a  semi-infinite  or  finite  edge  with  the  use  of  a  suitable  transition 
function  (in  general  transition  dyadic)  has  been  suggested.  Such  a  transition  func¬ 
tion  can  be  empirically  derived  by  requiring  the  continuity  of  the  total  field  at  the 
shadow  boundaries  of  the  direct  wave  (edge  wave),  namely  at  the  planes  4>i,2  = 
This  coj?tinuity  was  guaranteed  by  the  UTD  evaluation  of  edge  diffracted  field  as 
expressed  by  eqs.  (4.58a)-(4.58c),  but  it  is  violated  in  the  paraxial  region  after 
the  introduction  of  the  function  which  assures  the  uniformity  of  the  total 

diffracted  field.  Obviously,  the  edge  wave  can  be  multiplied  by  a  similar  transition 
function  so  that  the  total  field  retains  its  continuity  in  the  paraxial  region  as  well 
as  outside  of  it,  where  the  transition  function  receeds  to  unity.  In  addition,  such  a 
multiplicative  correction  would  yield  a  finite  total  field  along  the  extension  of  the 
edge. 

The  edge  wave  contribution  to  the  total  radiated  field  by  the  dipole- trihedron 
configuration  is  given  explicitly  by  eq.  (2.19)  which,  for  the  half  plane  case,  reduces 
to  eq.  (2.20).  For  the  plane  angular  sector,  a  convenient  modification  of  the  edge 
wave  reads: 

•|/3o  co8/JoMn(^/2)  +  ^  co8(^/2)] 
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(4.65) 


F^^I2]^k{iLt6\T?  gxp(-jA;o3o) 

y/sm0Q  so 

with  sq  denoting  the  distance  of  the  observation  point  from  the  dipole  and  is 

the  elevation  angle  of  the  observer  in  the  guiding  edge  fixed  coordinate  system 

centered  at  the  point  of  the  projection  of  the  dipole  onto  the  edge.  The  continuity 

of  the  field  at  the  boundary  ^  =  0,  is  shown  trivially  by  noting  that,  for 

(^  =  0,  /So  =  TT  -  a  - 

4.5  Discussion  and  numerical  results. 

The  edge  wave  vertex  diffracted  field  is  a  higher  order  term  with  respect  to 
the  large  parameter  in  the  asymptotic  solution  of  the  radiation  of  a  dipole  in  the 
vicinity  of  the  edge  of  a  trihedron.  Nonetheless,  it  contributes  significantly  to  the 
field  especially  in  the  paraxial  region  of  the  guiding  edge  and  along  its  extension. 
And  its  major  contribution  concerns  the  /S-directed  component  of  the  field.  The 
latter,  as  pointed  out  in  the  discussion  of  the  Physical  Optics  solution,  is  due  to 
the  accumulation  of  electric  current  flow  lines  in  the  vicinity  of  the  guiding  edge 
excited  by  the  ray  optical  edge  wave,  which  is  consistent  with  the  theoretically 
predicted  singularity  of  the  field  as  the  edge  is  approached. 

The  solution  based  on  Michaeli’s  equivalent  currents  is  essentially  an  asymp¬ 
totic  PTD  approach,  in  that  a  fringe  current  effect  due  to  the  terminating  edges  is 
added  to  the  uniform  edge  wave  currents.  The  rigorness  of  the  approach  may  be 
questioned  at  this  point,  since  the  derivation  of  the  fringe  edge  currents  assume  an 
infinite  edge  and  uniform  plane  wave  illumination.  However,  the  field  is  expected  to 
retain  its  singular  behavior  in  the  vicinity  of  the  vertex,  which,  moreover,  does  not 
contradict  the  “tip  condition”  (i.e.,  the  behavior  of  the  field  in  the  neighborhood 
of  a  vertex)  as  investigated  rigorously  by  several  authors.  Besides,  the  addition  of 
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the  fringe  current  to  the  physical  optics  current  effect  improves  the  total  solution 
as  shown  in  Figs.  29-34,  where  these  results  are  compared  with  moment  method 
ones.  It  should  be  emphasized  that  the  solution  attempted  in  this  chapter  by  no 
means  presents  a  complete  rigorous  representation  of  the  tip  diffracted  field,  but 
it  merely  includes  the  information  of  the  truncation  of  known  components  of  the 
currents  flowing  over  the  edge  surface  and  can  serve  as  a  good  engineering  approx¬ 
imation  to  the  problem.  Again  in  a  PTD  context  the  completeness  of  a  solution 
requires  the  addition  of  a  vertex  current  component,  i.e.,  a  current  excited  by  the 
vertex  of  the  trihedron,  which  however  remains  unknown  and  it  does  not  appear 
possible  that  it  can  be  extracted  from  the  rigorous  solutions  [1,2]. 

Comparisons  of  the  calcualted  field  (denoted  as  MEC  on  the  graphs)  with 
moment  method  results  (MM)  are  made  in  Figs.  29-34,  for  the  square  plate  shown 
in  Fig.  12.  The  direct  edge  wave  (EW)  properly  modified  in  accordance  with  the 
discussion  in  Section  4.4.3  is  also  plotted  so  that  the  effect  of  the  corner  diffracted 
fields  is  better  illustrated.  Note  that  in  Fig.  29  the  edge  wave  edge  diffracted  wave 
has  been  added  to  the  direct  dipole  field  (EW-J-EWED). 

The  field  evaluated  with  the  equivalent  current  method  shows  good  agree¬ 
ment  for  those  flat  plate  geometries  where  the  dipole  is  placed  in  the  close  vicinity 
of  one  of  the  edges  of  the  plate  and  sufficiently  far  from  its  corners,  especially 
in  regions  where  contributions  from  other  diffraction  mechanisms  other  than  the 
corner  adjacent  to  the  guiding  edge  are  known  to  be  negligible.  This  becomes 
evident  in  Figs.  29-34,  where  the  aRreemeiit  between  the  moment  method  data 
and  the  calculated  field  progressively  improves  for  larger  values  of  the  <f>  angle.  In 
fact,  for  <f>  =  60®,  other  mechanisms  such  as  diffraction  from  the  remote  corners 
of  the  square  plate  as  well  as  double  and  triple  edge  diffraction  may  contribute 
significantly,  while  their  effect  diminishes  as  »  180®.  Our  approximate  analysis 
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also  improves  when  the  paraxia]  region  is  approached  {0  — ►  0),  where,  as  a  matter 
of  fact,  the  total  field  is  stronger.  The  latter  justifies  the  validity  of  the  asymp¬ 
totic  analysis  which  resulted  into  the  multiplicative  introduction  of  two  transition 
functions  as  well  as  the  choice  of  the  correction  factor  for  the  direct  dipole  field. 
Unfortunately,  the  ^-directed  field  does  not  exhibit  an  analogous  agreement  basi¬ 
cally  because  the  tinalysis  does  not  include  secondary  mechanisms  which  clearly 
contribute  significantly  to  the  pattern  of  that  >  olarization. 

The  results  also  reveal  a  smsJl  variatio;r  of  the  total  field  with  respect  to  the 
azimuthal  coordinate  (angle  ^),  in  contrast  with  the  relative  large  changes  of  the 
calculated  pattern  in  the  elevation  plane.  The  4>  dependence  becomes  significant 
only  at  lower  cuts  where  the  contribution  from  the  opposite  edge  and  its  two 
adjacent  corners  is  appreciable. 

The  second  example  examined  involves  also  the  radiation  of  a  small  monopole 
in  the  close  vicinity  of  one  of  the  edges  of  a  rectangular  plate,  but  now  in  differ¬ 
ent  distances  from  its  corners.  The  geometry  of  this  monopole-rectangular  plate 
configuration  is  depicted  in  Fig.  35.  The  calculated  field  is  compared  again  with 
moment  method  results  as  well  as  measured  data,  as  shown  in  Figs.  36,  37  on  the 
azimuthal  planes  <t>  =  180®,  ^  =  150®,  respectively.  The  accuracy  of  the  measured 
data  deteriorated  for  measurements  in  smaller  azimuthal  angles,  where  the  support 
of  the  structure  influenced  significantly  the  measured  radiation  pattern.  In  Fig.  36, 
a  ram  has  been  placed  around  the  remote  corners  and  the  opposite  edge  of  the 
rectangular  plate,  so  that  their  effect  in  the  total  pattern  is  reduced.  Clearly,  the 
agreement  is  better  in  this  case  (Fig.  36),  in  contrast  with  the  results  of  Fig.  37 
where  the  ram  was  removed. 

It  should  be  noted  that  neither  the  Physical  Optics  solution  nor  the  equivalent 
current  formulation  are  expected  to  yield  accurate  results  for  small  angular  sector 
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angles  (a  <<  7r/2).  In  this  case  a  strong  coupling  between  the  two  edges  forming 
the  sector  occurs,  which  is  not  encountered  in  the  evaluation  of  the  radiation 
integral  of  the  currents  flowing  along  the  edge  (1).  Moreover,  when  a  — ♦  tt,  the 
Physical  Optics  corner  diffracted  field  vanishes  and  the  total  solution  reduces  to 
the  edge  wave  over  an  infinite  wedge,  whereas  the  same  property  is  not  true  for  the 
equivalent  current  formulation  of  the  vertex  diffracted  field.  The  latter,  therefore, 
fjiils  in  cases  of  very  wide  angles,  which  require  a  more  careful  treatment. 
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Figure  29:  Far  field  ^-directed  field  at  an  azimuthal  angle  ^  ss  fiO^  for  the  square 

plate  of  Fig.  12. 
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Fisure  30:  Far  field  /9-directed  field  at  an  asimutbal  angle  ^  s  90**  for  the  square 

plate  of  Fig.  12. 


70 


0.0  30.0  60.0  90.0  120.0  60.0  60.0 

ANGLE  (DEG.) 


Figure  31:  Far  field  /9-directed  field  at  an  azimuthal  angle  ^  s  120^  for  the 

square  plate  of  Fig.  12. 
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Figure  32:  Far  field  field  at  an  asimutbal  angle  ^  s  135^  for  the 

square  plate  of  Fig.  12. 
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Figure  34:  Far  field  $~^Tecitd  field  at  an  azimuthal  angle  ^  * 

zquare  plate  of  Fig.  12. 
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Figure  35:  Geonietrv  of  the  monopole-rectangular  plate  configuration  used  for 

comparison  with  measured  data. 
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Figure  37:  Far  field  /3-directed  field  at  an  azimuthal  angle  =  150®  for  the 
rectangular  plate  of  Fig.  35.  Solid  line:  Calculated  field,  Pwhcd  llflCi  Measur 
field,  Dotted  line:  Moment  method  solution. 


CHAPTER  V 


Conclusions 

The  major  objective  and  motivation  of  this  study  was  to  describe  approxi¬ 
mately  the  edge  wave  diffraction  mechanisms  associated  with  the  interaction  of  an 
edge  wave  and  the  vertex  of  an  ideal  trihedron  as  well  as  two  adjacent  vertices  in 
a  realistic  three  dimensional  structure  containing  finite  edges.  These  mechanisms 
may  contribute  significantly  in  the  radiation  patterns  of  several  antennas  (such  as 
rectangular  horns)  and  the  RCS  of  edge  structures  investigated  with  UTD  tech¬ 
niques,  and  therefore  their  incorporation  improves  the  total  UTD  solution  and 
extend  UTD  as  a  powerful  high  frequency  theory. 

Explicit  expressions  of  the  edge  wave  were  presented  in  Chapter  II,  derived 
from  the  canonical  half  plane  and  infinite  wedge  problems.  The  investigation  of 
such  ideaJized  configurations  was  dictated  by  the  need  of  a  complete  definition  of 
the  edge  waves  and  a  better  understanding  of  its  behavior.'  It  was  found  that  the 
edge  wave,  which,  independently  of  its  excitation,  is  essentially  an  edge  guided  ray 
optical  field,  highly  dominates  the  paraxial  region  of  the  edge  due  to  its  singular 
behavior  there.  These  characteristics  imply  that  there  must  be  a  strong  coupling 
between  antennas  placed  close  to  the  edge  of  a  wedge  (which  has  been  already 
proven  by  Buyukdura  [7])  and  between  two  adjacent  vertices  in  a  polyhedral  body, 
as  found  by  Sikta  [4]. 

The  next  simple  configuration  examined  was  the  trihedron  (a  truncated 
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wedge),  as  a  first  step  to  the  approximation  of  more  complex  geometries.  A  dipole 
radiating  in  the  close  vicinity  of  one  of  its  edges  produces  an  edge  wave  which 
behaves  essentially  as  a  ray  optical  field  in  the  vicinity  of  the  vertex,  provided 
that  the  distance  of  the  excitation  dipole  from  the  tip  is  sufficiently  large.  To 
approximate  the  vertex  diffracted  wave  a  Physical  Optics  approach  was  adopted. 
The  latter  was  improved  in  Chapter  IV,  by  adding  to  the  uniform  current  a  fringe 
current  component  flowing  in  the  vicinity  of  the  terminating  edges.  Although, 
the  approach  is  neither  rigorous  nor  complete  from  a  PTD  point  of  view,  it  yields 
comparable  results  with  moment  method  as  well  measured  data  and  can  be  used 
as  a  first  engineering  approximation  to  the  edge  wave  edge  and  vertex  diffraction 
problems. 

Within  the  limits  of  our  ray  optical  approximations  the  results  established 
in  Chapter  IV  can  be  possibly  utilized  for  the  examination  of  the  double  corner 
and  edge  edge  wave  diffraction  mechanisms,  with  several  practical  applications. 
The  latter,  being  also  a  possible  extension  of  the  present  UTD  so  that  it  may 
incorporate  higher  order  diffraction  mechanisms,  awaits  future  work. 
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APPENDIX  A 


A  Uniform  Asymptotic  Approximation  of  the  Integral 
4(*)  =  dt  for  »(p)  <  0 

In  the  above  integral  G{i)  is  a  slowly  varying  continuous  function  of  i.  More¬ 
over,  let  the  large  parameter  k  have  a  small  negative  imaginary  part,  i.e.,  it  is 
presumed  that  the  medium  of  propagation  is  slightly  lossy.  The  phase  function 
p(/)  is  stationary  at  the  point  ~  0).  Since  in  our  applications  g{t)  ex¬ 

presses  the  negative  sum  of  the  distance  from  a  point  Q  on  a  straight  line  in  space 
to  two  discrete  points  not  on  this  line  ,  with  t  denoting  the  signed  arc  length 
from  the  origin  (<  =  0)  to  the  point  Q,  then  g(t)  is  negative  definite  and  attains  a 
maximum  at  In  this  case  it  can  be  shown  that 

<7'(0)  =  |<7'(0)|  (A.l) 

and 

g''(t)  <  0  ;  for  every  real  t.  (A. 2) 

For  convenience,  it  is  assumed  that  eqs.  (A.l),  (A.2)  are  true  in  the  following 
analysis.  Besides,  the  final  result  can  be  easily  generalized  for  a  phase  function 
not  satisfying  the  above  imposed  properties. 

For  if  >  0  the  major  contribution  to  the  integral  Ip(k)  arises  from  the  vicini¬ 
ties  of  the  stationary  phase  (/  =  t,)  and  the  end  point  (t  =  0),  which  coincides 
with  the  branch  point  singularity  of  the  integrand.  We  examine  the  asymptotic 
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behavior  of  the  integral  Ip{k)  for  the  particular  case  where  p  =  — 1/2.  An  evident 
transformation  of  variables  is  then  (  =  y/it  which  yields  the  equivalent  form 


/-1/2(*')  =  2^  H{()txp{jkh{^)}di  (A.3) 

where  we  have  set 

=  M{)  =  »«’)  (A-4) 

Since  ^(t)  is  stationary  at  tj,  where  i,  is  real,  then  h{()  exhibits  three  collinear, 
equidistant  saddle  points  at  =  0*(l,2  =  (Note  that  (1^2  imaginary 

if  /(  <  0).  A  convenient  transformation  in  this  circumstance  is  [28] 

h{()  =  aa  -  (a  ^  (A.5) 

It  readily  follows  that 

ao  =  h{U,2)  =  9(t,)  (A.6) 


a  =  -sipn(<,)  2)  ~  M0)l  =  -s»pn(<,)  Iy/^(<»)  -  p(0)| 

Provided  that  a  is  not  very  large  the  integral  /_i/2(fc)  approximated  by 

/_l/2  2G{0)^\,=o  J^cxp{jk[a(i-{a  +  $^)^]}d3  (A.7) 

Using 

h"(0)  =  2p'(0)  ;  A"(ei)  =  4<,/(<,)  ;  h^*H0)  =  12 g'' {0)  (A.8) 

it  is  not  difficult  to  verify  that  the  mapping  derivative  at  s  =  0  for  a  sufficiently 
small  equals  to 
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\  is'(o)r  \  iw'(<,)r  [ip«(o)| 


(A.9) 


For  the  details  .of  the  derivation  of  eqs.  (A.8),  (A.9),  which  are  useful  in  the 
examination  of  the  limiting  behavior  of  the  asymptotic  solution  as  o  — »  0  (or 
-♦  0),  the  reader  is  referred  to  (28).  The  integral  involved  in  eq.  (A.7)  can  be 
expressed  in  terms  of  the  parabolic  cylinder  function  of  order  —1/2.  In  particular, 
one  obtains 

*xp{jfc(oo-tt^/2)}  X 


I>_l/2(«P(W4)av/2fc] 


We  can  employ  the  formula  (29) 


(A.IO) 


Pp(z)  =  exp(-ipjr)I>p(-r)  +  exp{-j(p  +  l)ir/2} I>«p_i(jr)  (A.ll) 

for  p  =  -1/2  and  rewrite  eq.  (A.IO)  as  follows: 

/-l/2(*f)-G(0) - - ^xp{jk(ao-a  /2)}^  • 


+  v/2exp(-j>/4)P_i/2(e3tp(-i»/4)lo|>/^]  ;<»  >0 


(A.12) 


1 /2  [«*P  O’t/^  )|a  I  v/^l 


;t$  <0 


It  is  reminded  that  the  above  representation  is  valid  for  small  values  of  the  param¬ 
eter  a.  An  approach  valid  for  larger  values  of  a  appears  in  general  cumbersome. 
However,  one  can  simply  perturbate  the  original  solution  given  by  eq.  (A.12),  to 
obtain  an  asymptotic  approximation  in  which  the  stationi^y  phase  contribution 
would  reduce  to  the  well  known  result  for  an  isolated  stationary  phase  point  m  a 
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becomes  larg«.  In  particular,  for  >  0  using  eq.  (A.9)  one  writes 


v/2exp(-jir/4)  -^j2^^2>_,/2lexp(-i7r/4)|a|v/2fc]}  (A.13) 

From  the  asymptotic  form  of  the  parabolic  cylinder  function  of  order  —1/2  (29], 
for  |a|v^I;  1,  one  obtadns 

Substituting  the  above  expression  into  the  second  term  of  eq.  (A.13)  it  yields  the 
isolated  stationary  phase  contribution 

V,  -  »/4]) 

y/U  ^k\g'\i,)\ 

A  transformation  similar  to  that  of  eq.  (A.5)  cannot  be  applied  in  the  case  of 
an  arbitrarily  negative  p,  because  the  resulting  integral  cannot  be  reduced  into 
closed  form  expression.  To  obtain  a  uniform  asymptotic  aproximation  we  adopt  a 
more  empirical  approach.  Namely,  it  appears  convenient  to  consider  a  quadrature 
approximation  of  the  phase  function  ^rft)  in  the  neighborhood  of  the  end  point  i  = 
0,  at  which,  as  mentioned  before,  the  integral  presents  a  branch  point  singularity 

9(0  ^  9(0)  +  9'(0)f  +  (A.14) 

Then,  based  on  this  approximation  the  behavior  of  the  end  point  effect  and  the 
coincident  branch  point  singularity  is  examined.  After  some  elementary  manipu¬ 
lation,  which  follows  the  same  steps  with  the  particular  case  examined  earlier  one 
finally  obtains  the  asymptotic  result: 

I^ik)  %  6’(0)r(-p)(2fc)P/2exp(ypir/4)  expOfc[9(0)  +  oV2]}  X 
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(A.15) 


where  now 


fjfp  (-ip»)  (i»/4)|o|v/2fc] 

^^p(exp(i»/4)|o|v^) 


\t,  >0 


;<#  <0 


v^2i<^”(oji 


(A.16) 


It  should  be  noted  that  the  above  equation  represents  the  end  point  effect  only. 
Evidently,  it  exhibits  an  abrupt  phase  change  of  exp(jp)r)  as  the  stationary  phase 
point  passes  through  the  origin.  This  discontinuity  should  normally  compensate 
the  discontinuity  of  the  stationary  phase  contribution  to  the  integral  Ip{k)  in  order 
for  the  total  asymptotic  solution  to  remain  uniform  across  the  boundsiry  tg  =  0. 
Such  a  stationary  phase  point  effect  which  guarantees  the  uniformity  of  the  total 
solution  is  derived  by  a  heuristic  multiplicative  correction  of  the  isolated  stationary 
phase  contribution  with  kiie  proper  transition  function. 

The  substitution 


=  exp(-jpir/4)(V^iair'’ exp(ifca^/2)I>p(exp(iir/4)|o|>/5ik] 

(A.17) 


furnishes  an  alternative  expression  for  the  end  point  contribution,  namely 

/®(t)  %  6'(0)r(-p)JfcPexp(-;>7r/2)exp{jM0)}[j'(0))Pi^(Ira*)  (A.18) 

where  the  branch 


(-1)^  =  exp(jpir) 

should  be  chosen.  In  the  above  equation  f^(’)  behaves  essentially  as  a  complex 
transition  function:  it  reduces  to  unity  for  large  values  of  its  argument  and  com¬ 
pensates  the  singularity  of  ^'(0)  as  the  stationary  point  approaches  the  end  point 
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at  f  =  0.  For  a  detailed  description  of  the  behavior  of  the  transition  function  Fp(*) 
the  reader  is  referred  to  Appendix  B. 


APPENDIX  B 


Edge  Wave  Transition  Functions 

Two  types  of  transition  functions  are  extensively  used  in  the  examination  of 
the  edge  and  vertex  edge  wave  diffraction  mechanisms  arising  from  the  integration 
over  the  singularity  of  the  edge  wave  which  coincides  with  the  end  point  of  the 
inlegration  path: 

•  The  edge  wave  vertex  transition  function 


•  The  edge  wave  edge  transition  function 

where  e  >  0.  The  latter  for  every  x  is  the  complex  conjugate  of  the  former,  which 
is  related  with  the  parabolic  cylinder  function  of  order  —v  via 

Fl^(x)  =  exp(ji'7r/4)(2|Tl)‘'/^exp(j|Tl/2)T)_^[exp(>ir/4)Y/2|x|]  (B.l) 

For  a  complete  discussion  of  the  properties  of  the  parabolic  cylinder  ^  •'.ction  the 
reader  is  referred  to  (29,30,31).  For  small  values  of  its  argument  it  assumes  the 


power  series  expansion 


(B.2) 


while  for  large  arguments  it  can  be  approximated  by  the  first  terms  of  the  asymp 


1) 


totic  expansion 

=  exp(-z2/4).*"'^ 
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Incorporating  expressions  (B.2),  (B.3)  in  the  definition  (B.l)  of  the  transition 
function  outside  the  transition  regions  one  may  assume  the  approximations 


«  >/»exp(ji/ir/4)l*r/*exp(>l*|/2) 
1 


r(*^)  r(|) 


2exp{iir/4)  ^  .  j(u  -  1/2)  . 

V  *  +  - _ *  + 


(B.4) 


for  |r|  <<  1,  and 


+  i/(i/  + l)(i/  +  2)(i/  +  3)  . 

Fi,{T)  ^  1  -  -  - --2^5 - + 


(B.5) 


Of  particular  importance  in  our  numerical  applications  is  the  transition  func¬ 
tion  of  order  -1/2  {v  =  1/2).  The  parabolic  cylinder  function  of  order  -1/2  can 
be  expressed  in  terms  of  the  modified  Bessel  function  K  of  order  1/4  [30,31]  as 
follows: 


=  (B.6) 

Incorporating  (B.6)  and  the  relationship  [30] 

into  eq.  (B.l )  for  i/  =  1/2  it  yields  the  expression 

=  5«<p(-j3»/8).xp(j|t|/2)yHij«j)>(|»l/2)  (B.7) 

Eq.  (B.7)  can  be  rewritten  in  terms  of  the  Bessel  functions  and  J_  namely 

=  ^exp(>»/4)«P(>kl/2) 

•  [exp(-iir/8)  J_,/4(|t|/2)  -  exp(j>/8)  Jj/4(1t1/2)] 

('b.8) 
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The  Bessel  functions  of  order  1/4  and  -1/4  are  tabulated  in  [32].  Therefore,  eq. 
(B.7)  can  be  readily  employed  for  the  numerical  computation  of  the  function 
xhe  relatively  simple  computer  program  that  follows  uses  a  quadrature  approxi- 

t 

mation  of  the  Bessel  functions  for  |t|  <  0.75,  an  approximation  up  to  order  |z|~^ 
for  |r|  >  2.5  and  a  third  order  Lagrangian  interpolation  in  the  transition  region,  to 
achieve  an  up  to  third  decimal  point  accuracy  of  the  edge  wave  transition  function. 
The  amplitude  and  the  argument  of  the  edge  wave  transition  function  is  plotted 
in  Figs.  38  and  39  respectively. 
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C  Subroutlnt  for  the  eveluetion  of  the  edge  wave 
C  transition  function  DCs). 

I 


COMPLEX. rUNCTION  pcrT(xo) 

I 

C  Nueber  of  points  in  the  Interpolation  region. 
PAXAHETER  INP-49 


I 


I 

C 

C 

c 

I 


c 

c 

c 


c 

c 


COMPLEX  CJ,CPll,CPd4,CPdl 
ZNTECEX  Z 

REAL  JM(lNP).JP(XNP),X(lNP).SM,SnLR.lC 
REAL  PI.C14.GJ4.jnl4,JPt4.«.XD.AM.lM,AP.BP 
REAL  DX1.0X2.P 

Define  Units  of  the  snail  argunent-linesr  lnterpolation>largs  argunent 
regions  for  the  spproxinstion  of  the  Bessel  Punctions  so  that  a  throe 
digit  accuracy  Is  achieved. 

XO-O.S*XO 

BHLR-l.B-S 

BH-.I 

BC-6. 

C14-3.fi25E0990l2 

C34-1.22S4167024 

PX-3.141S926S3 

CJ-(0.,1.) 

CPJ4.(. 707109711,. 70710C7B1) 

CPJS-( .9231979S12, .393613432) 

CPI9-CONJG(epJ9} 

DATA  %/ . 6,. B«l., 1.3, 3. 4,1. 6, X.B, 2. ,2.2,9. 4,2. 6,2. 9, 9. ,3. 2, 3. 4, 3. 6, 3. B 

1. 4..  4. 2. 4.4. 4. 6. 4. 6. 5..5. 2, $.4, 5.6,6. 6, 6. ,6. 2, 6. 4, 6. 6, 6.9,7., 7. 2 

2. 7. 4. 7. 6. 7. 6. 1. .  1. 2, B. 4, I. 6, B.i, 9. ,9. 2, 9. 4, 9. 6, 9. I, 10. ,10. 2/ 

DATA  OH/. 97369. .91704, .66931, .52599, .31619, .25116, .12290,  .00359, •.10447 

1, -.  19919, -.27977, -.34195, -.39751, -.41537, -.42557, -.41990, -.39622 

2,  -.35947, -.31059, -.25192, -.19609, -.11590, -.04397,. 02694,. 09403 

3..  155.. 20776.. 25064.. 26226.. 30166.. 309.. 10361.. 21665.. 25911.. 22195 

4..  17706.. 12639.. 07201.. 01615, -.039. -.091132, -.13993, -.17991. -.21291 
5, -.23671. -.25077, -.25466, -.24942, -.23251/ 

DATA  OP/. 75966,. 76901,. 75223,. 71291.. 65453,. 59035,. 49366,. 39791,. 29623 

1..  19233.. 09947. -.00916, -.10064, -.19235, -.2521. -.30915, -.34926 

2,  -.37476. -.36449, -.37965, -.35673, -.32551, -.26097, -.22722, -.16663 

3,  -.10174, -.03514,. 03057,. 09293,. 14967,. 19979,. 23966.. 269,. 296 

4..  29227.. 29699.. 27037.. 24363.. 20796.. 165.. 11653.. 06459.. 01129 
5, -.04131, -.09114, -.13633, -.17521, -.20639. -.22992/ 

Otep  site. 

B-.2 

Apprexinatien  of  the  Bessel  functions. 

XrCXD.lT.BH)  TBCN 
Basil  argunent  tern. 

XPIXD.LB.O.)  nCN 
PCPT»(0.,0.) 

BETURN 

BLBE  XP(XO.LT.BnLR)  TBCN 

PCfT-BQRTiH*B0»T(2.«XD))*CP09*CEXP(C0*XD)/C34 

BETURN 

SLSE 

3M1 4-60BT ( OORT ( 2 . /» ) ) • t 1 . -XD«XD/3 . eBD* • 4/4 2 . ) /634 
JPl 4-4 . •S0RT( f 0RT( . 5*XD) ) • ( 1 . -XD»XD/5 . eXD*  *4/90 . )/Gl 4 
BNDXr 

BLEB  Xr(XO.GE.Sn.AMD.XO.LT.BG)  TBEN 
interpolation  region. I,af range  three  point  interpolation  toraula  is 
utilised. 

BO  13  Z-2.1NP-1 

2r(XD.GB.B(Z).AIID.IO.LT.B(Xein  tRIN 
DXl-XD-S(Z) 
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DX2-XD-X(1^1) 

ir(ABS(DX2).CT.A»S(DXl))  OOTO  14 
ir(ABS(DX2).LE.A>S(DXl))  OOTO  IS 
XNoir 

13  CONTZNUK 

14  r-DXl/M 

jNi4-.S*p*(r-i. )*jn(z-i)«(i.-T*r)'*jN(X)«.S*r*(p«i.)*jn(iil) 

JP14-.5«P*(P-1.)«JP(I-1)^(1.-P«P)»JP(I)*.5‘P«(P*1.)*JP(X^1) 
GOTO  IB 

15  P»DX2/R 

jni4-.5*p*(p>i.)*jn(x)«(i.>p*p)*jn(x«i)4'.S*p*(p«l.)*JH(x«2) 
JP14-.5*P«tP-l.»*JPtX>-»(l.-p*p)*jp(l4l)*.5‘P»(P^l. J*JP(X#a» 
GOTO  IB 

C  LAi9t  •i9uBtnt  fora. 


ELSE  ir(XO.CE.»C)  THEN 
An-COS(XO-0.12S*PX) 

■H-SZN(XO-0.125*PX) 

A?-COS(XD-0.3?5»?!} 

BP>SIN(XD-0.375*PX) 

JM14-SQRT(2./(P1«XO))*(XH-»1.*M/(32.*XOI-105.*A«/(204B.»XD»XD)) 
aP14-SQXT(2./(TI*XOI )«(AP*3.*»P/(32.*XO)-105.*AP/(204B.«XD*XD) ) 
ENDXP 

IB  PCrT-SQAT(  PX«XO )  *CPJ4»  ( CUB* m4-.CPJB»  JP14  )  -CEXPC  CJ*XD) 

XETUXN 

END 
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Figure  38:  Amplitude  of  the  edge  wave  traniition  function. 
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